Mon. Not. R. Astron. Soc. OOP, (ipO] (2012) Printed 15 February 2012 (MN l^T^ style file v2.2) 



Berkeley Supernova la Program III: Spectra Near 
Maximum Brightness Improve the Accuracy of Derived 
Distances to Type la Supernovae 

^ Jeffrey M. Silverman, Mofian Ganesfialingam,^ Weidong Li/f 
O Alexei V. Filippenko 

^Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA 
^ Marc J. Staley Fellow 
^ ■[Deceased 12 December 2011 



o 

u 



Accepted . Received ; in original form 

ABSTRACT 

In this third paper in a series we compare spectral feature measurements to photomet- 
ric properties of 108 low-redshift {z < 0.1) Type la supernovae (SNe la) with optical 
spectra within 5 d of maximum brightness. The spectral data were obtained from 
Ph 1989 through the end of 2008 as part of the Berkeley SN la Program (BSNIP) and is 

Q presented in BSNIP I (Silverman et al., submitted) and the photometric data come 

}^ mainly from the Lick Observatory Supernova Search (LOSS) and are published by 

Ganeshalingam et al. (2010). The spectral measurements are presented and discussed 
Ci in BSNIP II (Silverman et al., submitted), and the light-curve fits and photometric pa- 

rameters can be found in Ganeshalingam et al. (in preparation). A variety of previously 
proposed correlations between spectral and photometric parameters are investigated 
^ using the large and self-consistent BSNIP dataset. We find the pseudo-equivalent width 

(pEW) of the Si II X4OOO line to be a good indicator of light-curve width, and the pEWs 
of the Mg II and Fe II complexes are relatively good proxies for SN colour. We also 
I employ a combination of light-curve parameters (specifically the SALT2 stretch and 

CNJ colour parameters xi and c, respectively) and spectral measurements to calculate dis- 

tance moduli. The residuals from these models are then compared to the standard 
f~>) model which uses only light-curve stretch and colour. Our investigations show that 

a distance model that uses Xi, c, and the velocity of the Si II X6355 feature does not 
^-H lead to a decrease in the Hubble residuals. We also find that distance models with flux 

ratios alone or in conjunction with light-curve information rarely perform better than 
the standard {xi,c) model. However, when adopting a distance model which combines 
the ratio of fluxes near ^3750 A and ^4550 A with both xi and c, the Hubble residuals 
^ are decreased by ^ 10 per cent, which is found to be significant at about the 2a level. 

The weighted root-mean-square of the residuals using this model is 0.130 ± 0.017 mag 
(as compared with 0.144 ± 0.019 mag when using the same sample with the standard 
model). This Hubble diagram fit has one of the smallest scatters ever published and at 
the highest significance ever seen in such a study. Finally, these results are discussed 
with regard to how they can improve the cosmological accuracy of future, large-scale 
SN la surveys. 

Key words: methods: data analysis - supernovae: general - cosmology:observations 
- distance scale 
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1 INTRODUCTION 

Type la supernovae (SNe la) have been used in the recent 
past to discover the accelerating expansion of the Universe 
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( Riess et al.|1998 Perlmutter et al.|1999 1, as well as to mea- 



sure cosmological parameters with increasing accuracy and 



precision (e.g., Astier et al.||2006 


iRiess et al.„2007| |Wood- 


Vasey et al.||2007 


Hicken et al. 


2009| 


Kessler et al.||2009| 


Amanullah et al. 


2010 


Suzuki 


st al. 


2011 1. In the most 



general terms, thermonuclear explosions of carbon/oxygen 
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(C/O) white dwarfs (WDs) are thought to give rise to SNe la 



(e.g., 


Hoyle & Fowler] 1960 Colgate & McKee 


1969 Nomoto 


et al. 


1984 see Hillebrandt & Niemeyer||2000 


for a review). 



However, after decades of observations and theoretical work, 
a detailed understanding of both the SN progenitors and ex- 
plosion mechanisms is still missing. In addition, there is very 
little known about how differences in the initial conditions 
in SNe la give rise to the measured range of observables. 
A large, self-consistent dataset is needed in order to solve 
these problems. 

The ability to do precision cosmology using SNe la re- 
quires that one is able to calibrate or standardise their lu- 



minosity. Phillips ( 1993 1 showed a tight correlation between 
light-curve decline rate and luminosity at peak brightness 
for the majority of SNe la, the so-called "Phillips relation." 
However, the addition of spectral observations to the light- 
curve data complicates the picture far beyond the simple 
assumption underlying the Phillips relation. Many compar- 
isons of spectral and photometric data of low-redshift SNe la 



have been performed in the past (e.g., Nugent ct al 



1995 



Benetti et al. 2005| Bongard e t al.| 2006[ Hachinger et al 



2006; Arsenij ivicet al.|2 008 Walker etal. 1201 l||Nordin et al. 
'2011b(|Blondin et al.pOllnChotard et al.|201ip . In addition, 
there has been similar work with SNe la at higher redshifts 



(e.g.,'Blondin et al."2006 



Altavilla et al.|2009[[Nordin et al' 



|201 lb. Walker ct al. ^ 2011 1. These studies were often aimed 
at finding a "second parameter" in SN la spectral or pho- 
tometric data which would increase the accuracy of their 
distance measurements. 

Most of these previous studies utilised relatively small 
and heterogeneous datasets. The data studied here, in con- 
trast, were self-consistently observed and reduced, and con- 
stitute one of the largest datasets to be analysed in this man- 
ner. Low-redshift {z ^ 0.1) optical SN la spectra from the 
Berkeley Supernova la Program (BSNIP; Silverman et al., 
submitted) are used along with complementary photometric 



data, largely from Ganeshalingam et al. (20101. The spec- 



tral features have been accurately and robustly measured 
(BSNIP II; Silverman et al., submitted), and the light curves 
have been fit using a variety of methods (Ganeshalingam 
et al., in preparation). 

We summarise both the spectral and photometric 
datasets used for this analysis in Section [2] and we describe 
our procedure for measuring spectral features, fitting light 
curves, and producing Hubble diagrams in Section [3] How 
these measured values correlate with each other and with 
previously determined classifications is presented in Sec- 
tion [4] along with our Hubble diagram results using various 
models for the distances to SNe la. We present our conclu- 
sions in Section [5] where the main results of our analysis 
are summarised and the most accurate and useful spectral 
indicators are discussed. Other, future BSNIP papers will 
expand on the analysis performed here with the addition of 
host-galaxy properties and late-time SN spectra. 



2 DATASET 

2.1 Spectral Data 

The same SN la spectral data are analysed in the current 
study as were used in BSNIP II. The spectra are all orig- 



inally published in BSNIP I. Most of the spectra were ob- 
tained using the Shane 3 m telescope at Lick Observatory 
with the Kast double spectrograph (Miller & Stone 1993j), 
and the typical wavelength coverage is 3300-10,400 A with 
resolutions of ~11 A and ~6 A on the red and blue sides 
(crossover wavelength ~5500 A), respectively. See BSNIP I 
For more information regarding the observations and data 
reduction. 

In BSNIP II, we required that a spectrum be within 
20 d (rest frame) of maximum brightness and we a priori 
ignored the extremely peculiar SN 2000cx (e.g. 
[2001), SN 2002 CX (e.g., [Li et al . 2003; Jha et al 



Li et al. 



SN 2005hk (e.g., |Chornock et al.|2006, Phillips et al 



and SN 2008ha (e.g., Foley et al.|2009 Valenti et al 



20061, 



20071, 



20091 



BSNIP II contains 432 spectra of 261 SNe la with a "good" 
fit for at least one spectral feature. However, only a subset of 
these data are used in the current study since not all of them 
have reliable photometric observations and we are currently 
only considering spectra within 5 d of maximum brightness. 
It was shown in BSNIP II that the spectral measurements 
do not evolve significantly during these epochs. For the 11 
SNe la that had more than one spectrum within 5 d of max- 
imum brightness and photometric information, we only use 
the spectrum closest to the date of maximum in the current 
analysis. 

For our investigation using arbitrary ratios of fluxes 



(see Section 4.101 we also use spectra within 5 d of max- 
imum brightness, even though previous studies only investi- 



gated spectra within 2.5 d of maximum ( Bailey et al.|2009 
Blondin et al.|[2bll l. When we employ spectra only in this 



narrower age range, our sample size for the flux-ratio analy- 
sis decreases from 62 to 38 objects. While our overall results 
are mostly unchanged when considering only spectra within 
2.5 d of maximum brightness, their signiflcance is weakened 
due to the smaller number of objects used. Adopting a larger 
age range than within 5 d of maximum increases the sample 
size only moderately and introduces much larger scatter in 
spectral measurements (see BSNIP II). 

2.2 Photometric Data 

A majority of the SNe in our spectral sample were discovered 
as part of the Lick Observatory Supernova Search (LOSS). 
LOSS is a transient survey utilising the 0.76-m Katzman 
Automatic Imaging Telescope (KAIT) at Lick Observatory 
( |Li et al.]|2000[ [Filippenko et al.||2001| see also Filippenko, 
Li, & Treffers, in preparation). KAIT is a robotic telescope 
that monitors a sample of ~ 15,000 galaxies in the nearby 
Universe (redshift z < 0.05) with the goal of finding tran- 
sients within days of explosion. Fields are imaged every 3- 
10 d and compared to archived template images, after which 
potential new transients are flagged. These images are exam- 
ined by human image checkers and the best candidates are 
re-observed the following night. Candidates that are present 
on two consecutive nights are reported to the community us- 
ing the International Astronomical Union Circulars (lAUCs) 
and the Central Bureau of Electronic Telegrams (CBETs). 
The statistical power of the LOSS sample is well demon- 
strated by the series of papers deriving nearby SN rates 
( Leaman et al.||2011[ |Li et al.||2011a|b I. 

In addition to the SN search, KAIT monitors active 
SNe of all types in broad-band BVRI filters. The first data 
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release of BVRI light curves for 165 SNe la along with de- 
tails about the reduction procedure have been published 



by Ganeshalingam et al. (20101. In summary, point-spread 
function (PSF) fitting photometry is performed on images 
from which the host galaxy has been subtracted using tem- 
plates obtained > 1 yr after explosion. Photometry is trans- 



formed to the Landolt system (Landolt 1983 19921 using 



averaged colour terms determined over many photometric 
nights. Calibrations for each SN field are obtained on pho- 
tometric nights with an average of 5 calibrations per field. 

We also include SN la light curves obtained from the 
literature to maximise the overlap between our photometric 
and spectroscopic samples. We include 29 objects from the 



Calan-Tololo sample (Hamuy et al. 19961, 22 objects from 
the CfAl sample ( |Riess et al.||1999| ), 44 objec ts from CfA2 
I Jha et ar]|2006 1 , and 185 objects from CfA3 (Hicken et al. 
2009 1 . In instances where we have data for the same SN from 



multiple samples, we use the light curve that is most densely 
sampled and best captures the light-curve evolution. We also 



include hght curves for SNe 1999 aw ([Strolger et al.||2002 |, 
1999ee, 2000bh, 2000ca, 2001ba (IKrisciunas et al.||2004al. 



2001bt, 2001cn, and 2001cz ( [Krisciunas et al.||2004b| ). Our 
final photometry sample consists of 335 multi-colour light 
curves, though we note that not all the objects in this sample 
have corresponding spectroscopy. 

Of the data within 5 d of maximum brightness inves- 
tigated in BSNIP II, 115 SNe have light-curve width or 
colour information and are included in the present study. 
The redshift range spanned by this sample is < z < 0.1. 
A complete list of these SNe la, their ages, spectral classifi- 
cations, and spectral feature measurements can be found in 
BSNIP II. The photometric parameters of these objects are 
presented by Ganeshalingam et al. (in preparation). 



3 MEASUREMENT PROCEDURES 
3.1 Spectral Measurements 

The algorithm used to measure each of nine spectral fea- 
tures and the features themselves are described in detail in 
BSNIP II. Here we give a brief summary of the procedure. 

Each spectrum has its host-galaxy recession velocity re- 
moved and is corrected for Galactic reddening (according 
to the values presented in Table 1 of BSNIP I), and it is 



smoothed using a Savitzky-Golay smoothing filter ( Savitzky 
|fc Golay||1964[ ). If the signal-to- noise ratio (S/N) is larger 
than 6.5 pixel"'^ over the entire spectral range, we attempt 
to define a pseudo-continuum for each spectral feature. This 
is done by determining where the local slope changes sign 
on either side of the feature's minimum. Quadratic func- 
tions are fit to each of these endpoints and the peaks of 
the parabolas (assuming that they are both concave down- 
ward) are used as the endpoints of the feature; they are then 
connected with a line to define the pseudo-continuum. We 
record the flux at the blue and red endpoints of the feature 
[Fb and Fr, respectively) as well as the pseudo-equivalent 
width (pEW; e.g., [Garavini et al.||2007[ ). 

Once a pseudo-continuum is calculated, a cubic spline 
is fit to the smoothed data between the endpoints of the 



spectral feature]^ From the wavelength at which the spline 
fit reaches its minimum (Amin) the expansion velocity [v) 
is calculated. The flux is then normalised to the pseudo- 
continuum, and the relative depth of the feature (a) and 
its full-width at half-maximum (FWHM) are computed. Fi- 
nally, every spectral feature in each spectrum is visually in- 
spected by more than one person and removed from the 
study if the spline fit and/or pseudo-continuum do not ac- 
curately refiect the spectral feature. 



3.2 Light-Curve Fitting 

A variety of methods have been developed to measure the 
photometric properties of SN la light curves. Here, we de- 
scribe three different light-curve fitting methods adopted in 
this paper to characterise the SN la light curves. While the 
light-curve parameters derived from each of these methods 
are degenerate to some extent, it is useful to perform all 
three fitting techniques for the purpose of comparing our 
results to previous results in the literature. There are also 
cases, for certain spectroscopic subtypes, where one method 
is superior to the other two. 



3.2.1 Template and Polynomial Fitting 

Our most direct measurement of light-curve properties 
makes use of the template-fitting routine introduced by |Pri-| 
eto et al. (20061. For a given photometric bandpass, a set 



of template light curves is used to construct models which 
match the light-curve data. The model light curves are linear 
combinations of the tem plate light curves using the weight- 
ing scheme described by Ganeshalingam et al. (20101. A x^- 



minimisation fitting routine is used to determine the combi- 
nation of templates that best fits the data. For band X, we 
measure the date of maximum brightness, the apparent peak 
magnitude (mx ) , and the light-curve width parametrised as 
the difference in magnitudes between maximum and fifteen 
days past maximum, Ami5(X). We independently fit the 
B- and V-band light curves for the SNe in our sample. 

In instances where we have a well-sampled light curve, 
but cannot achieve an acceptable fit with our template- 
fitting routine, we fit the data with a fourth-order polyno- 
mial. 

For both the template- and polynomial-fitting routines, 
we use a Monte Carlo routine to measure the uncertainty 
in our derived parameters. We simulate realisations of our 
dataset by randomly perturbing each data point by its la 
photometric error assuming a Gaussian distribution centred 
at mag. We fit the dataset realisation with our fitting rou- 
tine and measure light-curve properties for that simulation. 
We estimate the uncertainty in our derived parameters to 
be the standard deviation of 50 dataset realisations. 



^ No attempt is made to fit any of the Mg II or Fe II features in 
this manner; these complexes consist of so many blended spectral 
lines that it is unclear which reference wavelength to use when 
attempting to define an expansion velocity. 
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3.2.2 MLCS2k2 

The Multi-colour Light Curve Shape (MLCS) distance- 
fitting software was first introduced by Riess et al. ( 1996 1 



bundgut et al. 19931, since SALT2 was not developed to 



to simultaneously fit all light-curve data for a given SN la 
to produce a distance estimate. This method relies on the 
observation that more luminous SNe la have broader light 
curves ( Phillips|1993 1 and also have bluer colours during the 
photospheric phase. MLCS parametrises light-curve width 
using the parameter A, which measures the difference in 
absolute magnitude of the SN with reference to a fiducial 
SN la. MLCS attempts to disentangle intrinsic colour vari- 
ations from host-galaxy effects to also produce an estimate 
for host-galaxy extinction, Ay. MLCS2k2.v006 (referred to 
as simply MLCS2k2 for the rest of this work;|Jha et al. 12007 1 



is the the most current publicly available implementation of 
this fitting routine. In comparison to the original version, it 
has an expanded set of training templates and improvements 
in the treatment of host reddening and 7f -corrections. 

For our analysis with MLCS2k2, we use the Galactic 
line of sight prior which models the distribution of host- 
galaxy extinction values as a decaying exponential with a 
peak value of mag ( [Hatano et al.| |T 998| . We also set the 
host-galaxy Ry = 1.7 based on the cosmological analysis 



of Hicken et al. (20091. Their analysis found that a lower 



host-galaxy Ry reduced the scatter in the Hubble diagram 
compared to a more typical Galactic value of Rv ~ 3.1. 
A fit using MLCS2k2 is considered reliable (and thus its 
parameters are used in the current analysis) only when the 
reduced ^ 1-6. 



3.2.3 SALT2 

Spectral Adaptive Light-curve Template (SALT) was first 



developed by the SuperNova Legacy Survey (SNLS) (Guy 
et al.|[2005 |. Calculating distances with SALT is a two-step 
process. SALT first measures light-curve parameters that are 
expected to correlate with the intrinsic brightness of individ- 
ual SNe (i.e., light-curve width and colour). Then a model 
for the corrected apparent magnitude of a SN, ms,corr, is 
adopted which applies linear corrections for the light-curve 
width and colour to the measured apparent magnitude, mg. 
Thus, the corrected apparent magnitude has the form 

JTiB.corr = TTiB -|- a X (light-curvc width) — /3 X (colour). (1) 

The constants a and /3 are found by minimising using dis- 
tance estimates from a large sample of SNe la compared to a 
cosmological model (see Section [3.3| for more details). Unlike 
MLCS2k2, SALT avoids making assumptions on host-galaxy 
extinction and the intrinsic reddening of SNe. 

S ALT2 is an updated version of SALT with an expanded 
training set of light-curve templates and is the version im- 
plemented here. SALT2 is trained on light curves and spec- 
tra from low-2 SNe compiled from the literature and high-z 
SNe from the first two years of the SNLS ( |Guy et al.|2007 |. 
SALT2 measures a parametrisation of the light-curve width 
(a:i), the SN colour (c), and the apparent B-band magnitude 
at maximum light (ms). 

In fitting our light curves, we exclude /-band data, 
which are not included in the SALT2 template set. We also 
exclude subluminous SNe la, often of the spectral subclass 
of SN 1991bg-like objects (e.g., |Filippenko et al.|1992b[[Ler 



fit this particular subtype. This is achieved by using only 
SNe la with -3 ^ a;i sC 2 (as in |Blondin et al.||2011| ). Fi- 
nally, the results of SALT2 fits are utilised here only when 
the reduced x'^ < 2. 



3.3 Hubble Diagrams 

In this section we present the methodology used to standard- 
ise SNe la for cosmological application. We use a model that 
applies linear corrections for light-curve width and colour. 
The width of a light curve correlates with the intrinsic lu- 
minosity in the sense that SNe with broader light curves 
are also more luminous. This correlation has been well es- 



tablished ( [Phillips | 1993 1. The colour parameter combines 
the effects of intrinsic colour variations and host-galaxy red- 
dening. We use the SALT2 parameters xi and c as the 
parametrisations of light-curve width and SN colour, respec- 
tively. We will also generalise this approach to allow for lin- 
ear corrections using spectroscopic parameters. 

The distance modulus for each SN can be estimated 
from its redshift by n{z) = 25 + Slogj^g [Dl (z)], where Dl 
is the luminosity distance expressed in units of Mpc. The 
distance modulus including linear corrections for light-curve 
width and colour can be expressed as 



/tSN 



M - 



-13c. 



(2) 



The variables Q, and M (the fiducial absolute magnitude 
of a SN la) are determined by using a custom version of 
cosf itter (A. Conley, 2011, private communication) based 
on the Minuit minimisation package ( James fc Roos|1975 |. 
The software minimises the function 



X 



MSN,. 



(3) 



where ^(z) is the distance modulus of the galaxy in the 
cosmic microwave background (CMB) rest-frame redshift z, 
cTni is the measurement error in light-curve properties ac- 
counting for covariances between measured parameters, o-poc 
is the uncertainty due to deviations from Hubble's law in- 
duced by gravitational interactions from neighbouring galax- 
ies, and (Tint is a constant intrinsic scatter added to each SN 
to achieve a reduced ~ 1- We adopt 300 km s~^ as the 
peculiar velocity for each SN. The intrinsic scatter, (Tint, can 
be considered as the uncertainty associated with a model 
that attempts to standardise SNe using the parameters xi 
and c. 

Only objects with Zhciio > 0.01 are used in the Hub- 
ble diagrams in order to avoid including SNe with motions 
dominated by peculiar velocities. We also adopt the same 



Hubble diagram colour cut as Blondin et al. (2011) who ex- 



clude objects with c > 0.50. Finally, as mentioned above, 
all analysis using SALT2 fits is restricted to objects with a 
reduced <2. 

In this work we consider nearby SNe (median zcmb ~ 
0.021) and are not attempting to find a best-fitting cos- 
mology. A goal of this study is to combine photometric 
and spectroscopic properties such that SNe la become more 
accurate standardisable candles. We also aim to quantify 
the amount of improvement when using a variety of ob- 
served measurements. To that end, we adopt the standard 
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ACDM cosmology with Sim = 0.27, = 0.73, and w = -1 
when calculating /^(z). 



3.3.1 Models for Predicting SN la Distances 

Here we describe the generalisation of our calculation of the 
distance modulus for each SN to allow for linear corrections 
using measured spectral parameters (such as velocity and 



pEW; see Sections 4.2 and 4.4 1 or the ratios of pEWs and 
fluxes (such as 5R(Si II) and TZ; see Sections |4.5| and |4.10[ ). 
We consider five models for predicting distances to SNe la 
using combinations of the SALT2 measured light-curve pa- 
rameters {xi and c) and spectral measurements: 



MSN ~ 


TUB 


- M + 75, 


(4) 


MSN — 


niB 


— M + axi + jS, 


(5) 


MSN — 


TUB 


- M - I3c + ^S, 


(6) 


MSN — 


TUB 


- M + axi - Pc + jS, 


(7) 


MSN — 


niB 


- M + axi - Pc. 


(8) 



Here, S represents any spectral measurement: v, pEW, 
K(Si II), TZ, etc. The last model included in our study was 
already mentioned above (Equation[2| and is the model usu- 
ally adopted in cosmological studies of SNe la using only 
light-curve parameters (e.g., Astier et al.||2006 [Kowalski 
|et al.|[2008l [Hicken et aL]|2009| [Amanullah et"aLl(2010[ ). In 
the following analysis the so-called (xi , c) model will be the 
one to which we compare the other cosmological models that 
include spectral information. 



3.3.2 Cross- Validation 

Ideally, with a sufficiently large sample, the predictive abil- 
ities of a model could be inferred by inspecting the disper- 
sion of the residuals. However, for samples of limited size, 
the dispersion of the residuals is prone to statistical fluc- 
tuations and may not accurately reflect the true predictive 
nature of the model. Furthermore, for a flxed sample, one 
can always reduce the dispersion of the residuals by adding 
more variables to the model. However, it is not clear whether 
the added variables are actually improving the model itself 
or simply fltting to the noise inherent in the observables. 
For analysing the predictive nature of a model, it is useful 
to perform some form of cross-validation (CV) in which a 
subset of the entire sample is used to train the model and 
another subset is used to validate the predictive ability of 
that model. 



Bailey et al. ( 2009 1 use a sample of 58 SNe, 28 of which 
are adopted as a training set to train a model and the other 
30 are used as a validation set to assess the predictivity of 
the model. Using a smaller sample of 26 SNe, Blo ndin et al.| 
(20111) use a iC-fold CV method that allows all of the SNe 



to be used in the training and validation procedure. For this 



study, and following Blondin et al. (2011 1, we adopt the in- 
fold CV with K = 10. We also tested CV with K =2, 5, 
and 62 and find that our final results are mostly unchanged 
(however, see Section 4.10.21 for more on our K = 2 run, 
which is effectively what is used by Bailey et al.|2009 \ . The 
basics of the procedure are best illuminated by an example, 
as follows. 

Let us begin with a sample of 60 SNe. 10-fold CV starts 



by randomly dividing the sample into 10 subgroups of 6 SNe 
each. The first subgroup is set aside; this will be our first val- 
idation set. We combine the remaining 9 subgroups and train 
our model on the 54 other SNe to determine the best-fitting 
parameters (i.e., a, [3, 7, and M). Using the best-fitting pa- 
rameters found with the training set, we apply our model to 
the validation set and calculate the Hubble residual for each 
of the 6 SNe in our validation set. We repeat this process us- 
ing the second subgroup as the validation set and the union 
of the other 9 subgroups as the training set. This process 
is repeated a total of 10 times (once for each subgroup as a 
validation set) until we have calculated a residual for every 
SN in the sample. 

3.3.3 Comparing the Models 



As in Blondin et al. (2011 1, the dispersion in each model is 



estimated using the weighted root-mean-square (WRMS), 



WRMS^ = E.^i^4m(^.)-Msn,. 

where the weights, Wa, are given by 

_ 2 . 2 1 2 

The variance in WRMS is estimated as 

JV 



Var[WRMS] = [X]^ 



(9) 



(10) 



(11) 



The 1(7 uncertainty is found by taking the square root of 
the variance. This is a more appropriate estimator of the 
dispersion in the model then simply taking the standard 
deviation in the residuals since we are not guaranteed that 
the mean residual will be zero ([Blondin et al.|2011| . 

Following Blondin et al.l (2011), for each model which 



uses a spectral measurement (Equations |4||7| we also calcu- 
late the intrinsic prediction error (o-prcd), the intrinsic corre- 
lation {pxi.c) of the residuals with residuals using the (xi, c) 
model (Equation [8|, and the difference (Ajj^^c) in intrinsic 
prediction error with respect to the {x\, c) model. An uncer- 
tainty can be computed for ^xi,c (Appendix B of Blondin 
|et al.|pbll[ ), and thus the significance of the difference be- 
tween a given model and the standard (xi, c) model can also 
be computed. This parameter is the most direct comparison 
of how much better (or worse) a model which utilises a spec- 
tral measurement is compared to the {x\ , c) model, and how 
signiflcant the change is. Note that Axi,c < represents an 
improvement over the {xi,c) model. 



4 ANALYSIS 

4.1 Velocity Gradients 



Benetti et al. (20051 defined the velocity gradient, v = 
—Av/At, as the "average daily rate of decrease of the ex- 
pansion velocity" of the Si II A6355 feature and used this 
parameter to place each of their 26 SNe la into one of three 
categories. The high velocity gradient (HVG) group had the 
largest velocity gradients [v > 70 km d~^) and the low 
velocity gradient (LVG) group had the smallest velocity gra- 
dients. The third subclass (FAINT) had the lowest expan- 
sion velocities, yet moderately large velocity gradients, and 
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consisted of subluminous SNe la with the narrowest hght 
curves (Ami5(_B) > 1.6 mag). 

Even though the BSNIP sample is not well suited to ve- 
locity gradient measurements (the average number of spec- 
tra per object is ~2, shown in BSNIP I), we are still able to 
calculate a i; value for many of our SNe la. Figure [l] shows 
44 objects and their i) measurements plotted against their 
Ami5(_B) values. The points are colour-coded by their near- 
maximum Si II A6355 velocity, with red points representing 
normal-velocity objects and blue points representing high- 
velocity (HV) objects. These so-called "Wang types" were 
first presented by |Wang et al.| (|2009|, and in BSNIP II we 
discuss our definition of these subclasses in more detail. The 
data are also shape-coded by the aforementioned "Benetti 
types" (FAINT are stars, LVG are squares, and HVG are 
triangles). The top panel of Figure [l] shows all objects for 
which a « and Ami5(_B) are measured, while the bottom 
panel shows a close-up view of the same data such that the 
axis ranges match those of Figure 3b of [Benetti et al.| ( |2005[ ) . 

As a result of the definitions of the different subclasses, 
the FAINT SNe are found to the right in Figure[l](i.e., large 
Ami5{B) values) and the HVG SNe are found in the upper 
part of the Figure (i.e., large v values), though there is no 
obvious break between the various classes. As pointed out 
in BSNIP II and confirmed in Figure [ij the HVG and LVG 
objects have similar average Amii,{B) values and ranges of 
values. 



It was stated by Benetti et al. ( 2005 1 that v is weakly 



correlated with Amis{B), though there is no evidence of 
such a correlation in the BSNIP data. They also claim that 
there are three distinct families of SNe la (LVG, HVG, and 
FAINT) based partially on their plot of w versus Ami5(S). 
With almost 70 per cent more objects, the BSNIP data fill in 
this parameter space and cast serious doubt on the existence 
of truly distinct families of SNe la based on velocity-gradient 
measurements (see Table 4 of BSNIP II for the median val- 
ues of V and Ami5{B) for each of these three subclasses). 



4.2 Expansion Velocities 

4-2.1 Velocities at Maximum Brightness 

Expansion velocities of SNe la are calculated from the min- 
ima of various absorption features (see BSNIP II for more 
information on how this measurement is performed on the 
data presented here). These velocities have been compared 
to light-curve width measurements and photometric colours 
in a variety of ways (e.g. 



Hachinger et al. 2006 Blondin 



et al.|2011||Nordin et al.|2011b[ ). As discussed in BSNIP II, 



Hachinger et al. (20061 interpolate/extrapolate their ex- 



pansion velocities to the time of maximum brightness (i.e., 
t = d), and vq was defined as the expansion velocity of Si II 
A6355 at maximum brightness. They then compare these ve- 
locities to the light-curve shape parameter Ami5(_B). Fig- 
ure |2] presents the 44 SNe in the BSNIP data for which both 
vo and Ami5(_B) are calculated. As above, the points are 
colour-coded by "Wang type" and shape-coded by Benetti 
type. 

As in Figure [l] the FAINT objects are (by definition) 
found at the right edge of Figure[2] Similarly, the HV objects 
are all in the upper half of the Figure. All of the objects, ex- 
cept for a few of the highest velocity SNe and the object with 
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Figure 1. The velocity gradient versus Ami5(_B) of 44 SNe (top) 
and a close-up view of the low-i) objects (bottom). Blue points are 
high-velocity (HV) objects, red points are normal-velocity ob- 
jects, and black points are objects for which we could not deter- 
mine whether the SN was normal or high velocity (see BSNIP II 
for further details regarding how HV SNe are defined). Stars 
are FAINT objects, squares are low velocity gradient (LVG) ob- 
jects, and triangles are high velocity gradient (HVG) objects (see 
BSNIP 11 for further details regarding how these subclasses are 
defined) . 



the lowest velocity in the Figure, are within ~1500 km s~^ 
of 11,000 km s~^. This is remarkably similar to what was 
found by iHachinger et al. 



(20061, except the BSNIP data 



have a slightly la rger scatter around the typical velocity of 
~11,000 km s~\ [Hachinger et al.| ( |2006| ) also note that the 



majority of the scatter in velocity comes from the HVG ob- 
jects, which is also true in Figure |2] However, there is a ma- 
jor difference between the two results. The BNSIP data in 
Figure [2] show that HVG SNe have a huge range of vq, span- 
ning well above average to significantly below average val- 



ues, whereas the data presented by Hachinger et al. ( 2006 \ 



exhibit evidence of the (oft-quoted) one-to-one relationship 
between HVG and HV SNe la. As mentioned in BSNIP II, 
while most HVG objects are found to have expansion veloc- 
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Figure 2. The velocity at maximum brightness of the Si II A6355 
feature versus Am,i5(i?) of 44 SNe. Colours and shapes of data 
points are the same as in Figure [l] 



ities above the LVG objects, this is not an exclusive feature 
of HVG SNe la. 

A relationship between the calculated intrinsic "pseudo- 
colour" -Bmax — Knax (i-B., the _B-band magnitude at B-band 
maximum minus the l/-band magnitude at V'-band maxi- 
mum corrected for host-galaxy reddening) and has been 
seen previously (Foley et al. 2011 \. Figure |3] shows 39 BSNIP 
SNe for which we measure vo and the difference between ob- 
served B-band magnitude and observed \/-band magnitude 
at the time of B-band maximum brightness (what we refer 
to as {B — in this work). We opt to use (observed) 

[B — V)^^^ since it is an actual, physical colour of the SN 
at a discrete period of time, though keep in mind that in 
this study we do not attempt to correct for host-galaxy red- 



dening. Following Foley et al. (2011 1, we are only presenting 
SNe la with (B - F)max < 0.319 mag in Figure [s] 

The linear least-squares fit to all of the points is shown 
by the solid line. We see effectively no evidence for any over- 
all correlation (correlation coefficient of 0.17), and the lin- 
ear fit seems to be driven solely by the two outliers with 
(B— y)max > 0.25 mag. If those two objects are removed, the 
correlation coefficient drops to 0.05 (the best-fitting line to 
the remaining data is shown by the dotted line) . Finally, the 
dashed line is the relationship between vq and Bmax — Vmax 



from the model spectra of Kasen & Plewa (2007), as shown 
by Foley & Kascn (2011). When including all of the data 
in Figure [3| the correlation between vo and (B — V)^ 



weaker that what was seen by Foley et al. ( 2011[ ), where they 
derive correlation coefficients of 0.28 and 0.39 for two differ- 
ent datasetsj^However, no correlation is present whatsoever 
if the two significant outliers are removed. The difference be- 



tween the results shown here and those of Foley et al. ( 2011 1 



is likely due to the fact that we do not correct the BSNIP 
colours for any possible host-galaxy reddening while |Foley| 



The first dataset, yielding a correlation coefficient of 0.28, 



comes from |Wang et al.| (|2009| and includes some of the datal 
used here. However, the two samples are distinct and the method 
of velocity determination is significantly different. 



Figure 3. The velocity at maximum brightness of the Si II A6355 
feature versus observed (B — V)^^^ (i-e., the difference between 
observed _B-band magnitude and observed V-band magnitude at 
the time of _B-band maximum brightness) of 39 SNe. Colours and 
shapes of data points are the same as in Figure [l] The solid line 
is the fit to all of the data while the dotted line is the fit only 
to objects with (B — y)max < 0.25 mag. The dashed line is the 
relationship between vq and Bmax — Vmax from the model spectra 
of |Kasen fc Plewaj ( [2007| , as shown by [Foley fc Kasen| l [20lTl l . 



et al. (20111 attempt to convert the observed colours into 



intrinsic colours. We will delve deeper into this colour con- 
version in future BSNIP studies (Ganeshalingam et al., in 
preparation) . 



4-2.2 Velocities Near Maximum Brightness 

If instead of no we plot the actual measured velocity of the 
Si II A6355 feature for each object having a spectrum within 
5 d of maximum brightness versus Ami5(B), the same ba- 
sic trends are seen (but with nearly twice as many data 
points). A comparison of the velocity of the Si II A5972 fea- 
ture (within 5 d of maximum brightness) with ATni5(B) 
yields nearly identical results. The biggest difference is that 
the velocities are clustered around 10,300 km s~^, lower than 
that of the Si II A6355 feature. This difference between these 



features has been seen in previous studies as well ( Hachinger 
et al.|2006[ ). The same analysis using the velocity of the S II 
"W" once again shows the same behaviour, but with an even 
lower typical velocity (~9000 km s~^). This has also been 
pointed out in earlier work ( Hachinger et al.|2006 |. The ve- 
locity of the S II "W" feature is further discussed below. 

[Blondin et al.| ( |2011| ) present a Hubble diagram that 
is corrected by SALT2 light-curve width parameter xi and 
colour parameter c in addition to the velocity of the Si II 
A6355 feature. This yielded approximately a 10 per cent de- 
crease in the scatter of their Hubble diagram. It has also 
been shown that the Si II A6355 velocity is uncorrelated 
with both x\ and c, and thus it gives information beyond 
light-curve width and colour. However, the anticorrelation 
of this velocity with Hubble residuals (corrected for light- 
curve width and colour) is relatively small ( [Blondin et al.| 
20TT1). 
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Plotted in Figure |4] are the 66 SNe la in the BSNIP 
sample which have SALT2 fits and measured Si II A6355 
velocities within 5 d of maximum brightness. The velocities 
are plotted against xi, c, and Hubble residuals corrected 
for light-curve width and colour (only for objects which are 
used to make the Hubble diagram). As above, the points 
are colour-coded by "Wang type." In Figure [4] the data are 
shape-coded by what is referred to in BSNIP II as the "SNID 
type." The SuperNova IDentification code (SNID; Blondin 
& Tonry 20071, as implemented in BSNIP I, was used to 
determine the spectroscopic subtype of each SN used in 
BSNIP II (as well as this study). SNID compares an in- 
put spectrum to a library of spectral templates in order to 
determine the most likely spectroscopic subtype. Spectro- 
scopically normal objects are objects classified as "la-norm" 
by SNID. 

The spectroscopically peculiar SNID subtypes used here 
include the often underluminous SN 1991bg-like objects 
("Ia-91bg," e.g., Filippenko et al.|ri992b Leibundgut et al 



1993), and the often overluminous SN 1991T-hke objects 



("Ia-91T," e.g., Filippenko et al.|1992a 


Phillips et al. 


19921 


and SN 1999aa-like objects ("Ia-99aa," 


Li et al.|2001 


Strol- 



information regarding our implementation of SNID and the 
various spectroscopic subtype classifications. In Figure |4j 
la-norm objects are plotted as squares, Ia-91T objects are 
shown as upward-pointing triangles, Ia-99aa objects are dis- 
played as downward-pointing triangles, and objects with no 
determined SNID subtype are circles. 

The data in Figure |4] are extremely similar to those of 
Blondin et al. ( 2011 1 and show the same correlations (or lack 



thereof). The velocity of the Si II A6355 feature as measured 
from the BSNIP data is uncorrelated with xi (correlation 
coefficient of only —0.21). This velocity is also uncorrelated 
with c, as seen in Blondin et al. (2011 1. Removing SNe with 



c > 0.5 from the BSNIP data (which is done when Hubble 
diagrams are produced using this parameter) yields a cor- 
relation coefficient of 0.13 between Si II A6355 velocity and 
c. 

The bottom panel of Figure [4] shows the velocity of 
the Si II A6355 feature versus the {xi , c)-corrected Hubble 
residuals. If the correction term in a given model (here, the 
velocity of Si II A6355) is well correlated with the SALT2- 
corrected Hubble residuals, then the extra term is likely pro- 
viding new information that is actually in the data. Thus, 
the model is improving the fit not by fitting to noise, but 
to physical information contained in the data. However, the 
correction term here is uncorrelated with the residuals (cor- 
relation coefficient of 0.19). 

No improvement is found (i.e., the WRMS increases) 
when adding the Si II A6355 velocity to the standard {xi,c) 
model (A^i,c = 0.0083±0. 0084). [Blondin et al.| ( [20Tl| ) found 
that there was a < 10 per cent decrease in the WRMS when 
using the xi, c, and Si II A6355 velocity, but their Aa,^,c is 
consistent with (as is ours). They also find only a "modest" 
correlation between Si II A6355 velocity and (ii, c)-corrected 
residuals (correlation coefficient of 0.4; Blondin et al.|20"TT |. 
Therefore, it seems that adding the Si II A6355 velocity to 
the standard {xi, c) model does not significantly improve the 
precision of SN la distance calculations. 

In Figure [5] we plot all measured velocities of the Si II 
A6355 feature (within 5 d of maximum brightness) against 
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Figure 4. The velocity of the Si II A6355 feature versus SALT2 
light-curve width parameter xi (top), SALT2 colour parameter 
c (middle), and Hubble residuals corrected for light-curve width 
and colour (bottom). Blue points are HV objects, red points are 
normal- velocity objects, and black points are objects for which we 
could not determine whether the SN was normal or high veloc- 
ity. Squares are la-norm, upward-pointing triangles are Ia-91T, 
downward-pointing triangles are Ia-99aa, and circles are objects 
which do not have a SNID subtype (see BSNIP I for further de- 
tails regarding how these subclasses are defined). 
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[B — This differs from Figure [s] since the veloci- 

ties in that plot were interpolated/extrapolated to t = d 
(i.e., Wo) and this plot shows the actual velocities measured 
from the spectra. The top panel shows all 77 SNe from the 
BSNIP data for which both of these values have been mea- 
sured and the bottom panel shows a close-up of objects with 
[B — V)^„^ < 0.319 mag (in order to match the sample fit 



by Foley et al.|2011i 



The linear least-squares fit to all of the points is shown 
by the soUd line and the fit to SNe with {B — V)^^^ < 
0.319 mag is shown by the dotted line; the correlation co- 
efficients are 0.23 and 0.30, respectively, and slightly lower 
than the value found by |Foley et al. (20111, 0.39 (though 



the linear fit shown here matches well to what was found in 
their study). As with Figure [3j there is only marginal evi- 
dence for a correlation and there is a large amount of scatter 
around the linear fit. Again, this is unsurprising since we are 
measuring observed (B — V)^^^ and Foley et al. ( 2011[ ) are 
plotting pseudocolours that have been corrected for host- 
galaxy reddening. 

The typical {B — V)^^^ for HV objects is larger than 
for la-norm and Ia-91T/99aa, but smaller than that of la- 
91bg (which appear in the Figure as highly reddened, signif- 
icant outliers). However, the range of (B — V^),nax spanned 
by la-norm and HV objects is similar, with a significant 
amount of overlap. The scatter in {B — V)^^_^ of the HV 
objects (0.098 mag) is effectively equal to that of the other 
objects in the bottom panel of Figure [5] (0.101 mag). While 
this matches the scatter in the HV objects found previously 



(0.095 mag; [Foley et al.||2011[ ), it differs from the SNe with 
velocities < 11,800 km s~^, which they found to have a 
smaller intrinsic colour scatter (0.072 mag). The correction 
for host-galaxy reddening applied by Foley et al. (20111 is 
likely responsible for the decrease in colour scatter of the 
normal-velocity objects. 

The dashed line in Figure |5] is the relationship between 
vo and Bmax — Vmax from the mod el spectra of [Kasen fc| 



Plewa| ( |2007[ ) , as shown in Figu re 8 of [Foley fc Kasen] ( |2011 1 
Interestingly, even though ^Foley et al.| ( |2011[ ) plot intrinsic 
colours and we plot observed colours, both studies match 
these predictions extremely well. One difference between 
Figure [s] and Figure 8 of [Foley fc Kasen| ( [2011[ ) is that 
the BSNIP data contain a handful of objects that are ex- 
tremely reddened (i.e., they have relatively large values of 
{B — V)^^^). However, this can easily be explained. All of 
the la-norm and HV SNe in Figure [s] with {B - V)^^^ > 
0.31 mag have been observed to have significant reddening 
from their host galaxies (which is not taken into account in 
the models of ^Kasen fc Plewa|2007 K The other objects with 
{B — V)^^^ > 0.31 mag are Ia-91bg, which were also not 
Kasen fc Plewa| ^007\ 



discussed in the models of 



Figure 8 of Foley fc Kasen (2011 1 also presents the the- 



oretical relationship between velocity of the Ca II HfcK fea- 
ture and intrinsic {B — V)^^^. A comparison of this velocity 
and observed {B — V')^^;^,^ measured from the BSNIP data 
shows a very weak correlation with a correlation coefficient 
of only 0.14. Again, the scatter in {B — V)^^^ is similar for 
the HV and normal-velocity objects. However, this is unsur- 
prising since, as pointed out in BSNIP II, the Ca II HfcK 
velocities of HV and normal-velocity objects (determined 
using the Si II A6355 velocity) are highly overlapping. 

A distance model involving xi, c, and the velocity of the 
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Figure 5. The velocity of the Si II A6355 feature versus observed 
{B — V)„-jax ^^'^ ^ close-up of objects with {B — V)^^^ < 

0.319 mag {bottom). Blue points are HV objects, red points are 
normal- velocity objects, and black points are objects for which we 
could not determine whether the SN was normal or high velocity. 
Squares are la-norm, stars are Ia-91bg, upward-pointing triangles 
are Ia-91T, downward-pointing triangles are Ia-99aa, and circles 
are objects which do not have a SNID subtype (see BSNIP I for 
further details regarding how these subclasses are defined). The 
solid line is the fit to all of the data while the dotted line is the fit 
only to objects with [B — V)max < 0.319 mag. The dashed line 
is the relationship between vq and Bmax — Vmax from the model 
spectra of [Kasen fc Plewa[ l |2007| , as shown by [Foley fc Kasen] 
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Ca II HfcK feature was calculated, and while the WRMS 
technically decreased with the addition of this velocity, it 
was not found to be significant {A^j^^c ~ —0.0085 ± 0.0141). 
Other models involving this velocity were all found to de- 
grade the accuracy of distance measurements when com- 
pared to the standard (xi,c) model. 

Furthermore, we used Equations [4}|7] along with veloc- 
ities of all seven spectral features for which velocities were 
measured and compared the results to the {xi,c) model. 
The vast majority of these models predicted a larger scat- 
ter than the standard model corrected for light-curve width 
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and colour. However, both the O I A7773 triplet and the 
Ca II near-IR triplet, when combined with xi and c, were 
found to perform equally as well as when using just xi and 
c. Thus, adding either of these velocities did not degrade the 
distances calculated, but they did not significantly improve 
them either. 

On the other hand, the velocity of the S II "W," when 
used in conjunction with xi and c, decreased the WRMS by 
~3 per cent and the (Tprod by ~I4 per cent, at the 1.8a level 
(A^i,c = -0.0119 ± 0.0066). Figure |6] shows the 64 SNe la 
in the BSNIP sample which have SALT2 fits and measured 
S II "W" velocities within 5 d of maximum brightness. The 
velocities are again plotted against xi, c, and Hubble resid- 
uals corrected for light-curve width and colour (only for SNe 
that are part of the Hubble diagram). 

Neither xi nor c show any correlation with the velocity 
of the S II "W" (correlation coefficients of 0.13 and 0.17, 
respectively). Even when removing SNe with c > 0.5 from 
the BSNIP data (as done for the Hubble diagrams), the cor- 
relation coefficient becomes only —0.15. The bottom panel 
of Figure |6] shows the velocity of the S II "W" versus the 
(xi, c)-corrected Hubble residuals. The correction term is 
uncorrelated with the residuals (correlation coefficient of 
0.003). Figure [t] shows the actual Hubble diagram residu- 
als for this model (top panel) as well as the standard {xi, c) 
model (using the same set of objects; bottom panel) versus 
redshift. Also shown, as the grey band, is the WRMS for 
both models. 

While the relative depth of this feature has been seen 
to improve Hubble diagrams ( Blondin et al.|2011 and Sec- 



tion |4.3 of this work) , the velocity of this feature has not 
previously been shown to do this. When adding the velocity 
of the S II "W" feature to the standard {xi,c) model, the 
overall decrease in WRMS is relatively small, but the effect 
appears to be fairly significant. This distance model should 
be explored further using future, larger datasets. 



4.3 Relative Depths 

The depth of the bluer absorption of the S II "W" feature rel- 
ative to the pseudo-continuum has been shown to decrease 



the scatter of Hubble residuals by about 10 per cent ( Blondin 
|et al.|20TT| ). The relative depth (a) of this feature was found 
to be uncorrelated with both xi and c by [Blondin et al.| 
(20111, and its correlation with Hubble residuals (corrected 
for light-curve width and colour) is relatively small. Figurejs] 
presents the 64 BSNIP SNe la which have SALT2 fits and 
measured relative depths of the redder absorption of the S II 
"W" feature within 5 d of maximum brightness. The depths 
are plotted against xi, c, and Hubble residuals corrected for 
light-curve width and colour (for objects that are in the Hub- 
ble diagram). It should be noted that whereas [Blondin et al.| 
(2011 1 measure the bluer absorption of this feature (A5454), 
we measure only the redder absorption (A5624) in BSNIP II. 
While these absorptions are separated by < 200 A, there 
may be differences between the relative depths of the two. 
Nevertheless, we will compare the results presented here to 
those of [Blondin et al.| ( |2011[ ) with the caveat that we may 
be comparing "Red Delicious apples" to "Granny Smith ap- 
ples." 

The BSNIP data show a fairly significant correlation be- 
tween the relative depth of S II "W" and xi (a correlation 
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Figure 6. The velocity of the S II "W" feature versus SALT2 
light-curve width parameter xi, SALT2 colour parameter c, and 
Hubble residuals corrected for light-curve width and colour. 
Colours and shapes of data points are the same as in Figure [5] 
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Figure 7. Hubble diagram residuals versus zcMB for the 
(xijCjS II "W" velocity) model (top) and the standard (xi,c) 
model (Equation [s] bottom). The grey band is the WRMS for 
each model. Colours and shapes of data points are the same as in 
Figure [5] 



coefficient of —0.52). This is in stark contrast to I Blondin] 
et al. (20111, who found no evidence of such a correlation 
(correlation coefficient of —0.04). As opposed to xi, both 
studies agree that c is uncorrelated with the relative depth 
a of the S II "W" feature (the BSNIP data having a correla- 
tion coefficient of 0.09). The Hubble residuals corrected for 
xi and c show weak evidence of a correlation with the rel- 
ative depth of the S II "W" (correlation coefficient of 0.29). 

ondin et al. (20111 found that a model which includes 
xi, c, and a of the S II "W" will decrease the WRMS by 
< 10 per cent, while we find the WRMS to be effectively 
unchanged whether or not one adds in the relative depth of 
S II "W" (A^i,c = 0.0012 ± 0.0047). 

All other spectral features' relative depths were used, 
along with Equations |4||7| to create Hubble diagrams. No 
model significantly decreased the residuals over the stan- 
dard {x\,c) model. Models involving the relative depth of 
Ca II H&K, Si II A6355, and the O I triplet, each in combi- 
nation with xi and c, were found to be as accurate as the 
(xi, c) model. 
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Figure 8. The relative depth of the S II "W" feature versus 
SALT2 light-curve width parameter xi {top), SALT2 colour pa- 
rameter c (middle), and Hubble residuals corrected for light-curve 
width and colour [bottom). Colours and shapes of data points are 
the same as in Figure [s] 
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4.4 Pseudo-Equivalent Widths 



Nordin et al. ( 2011b I fit the temporal evolution of their pEW 
measurements to attempt to "remove" the age dependence 
of the pEW values. To do this, an epoch-independent quan- 
tity called the "pEW difference" (ApEW) was defined; it is 
simply the measured pEW minus the expected pEW at the 
same epoch using the linear or quadratic fit. In BSNIP II 
we calculated ApEW for the BSNIP sample. However, the 
relationships seen in BSNIP II involving ApEW values were 
also seen when simply using the pEW values (within 5 d of 
maximum brightness). This is due to the fact that pEWs 
do not evolve much within a few days of maximum. Fur- 
thermore, the ApEW values rely on defining a fit to the 
measurements which adds another assumption to the analy- 
sis. Thus, the current study will focus solely on pEW values 
within 5 d of maximum brightness and will not further in- 
vestigate ApEW values. Note that comparisons to |Nordin| 
et al. (2011b I will be made, despite the fact that their study 



uses ApEW values almost exclusively. 



4.4.1 Si 11X4000 



The pEW of the Si II A4000 feature has recently been found 
to be an indicator of light-curve width due its relatively tight 
anticorrelation with the SALT2 xi paramete r ([Arsenijevic 



et al. [1 2008: "Walke r et al.||2011| [Blohdin et al.||2011[ [Nordin 



et al.poilb, ^Chotardet al.||2011[ ). Curiously, in BSNIP II, 



only a weak correlation was found between this pEW and 
another often-used SN la luminosity indicator, the so-called 
"Si II ratio", 5R(Si II) (originally defined by Nugent et al. 
1995; see also BSNIP II and Section[43]for more information 
on this spectral parameter). Here the pEW of Si II A4000 is 
compared directly to photometric parameters. 

In Figure |9] we present the 57 BSNIP SNe which have 
SALT2 fits and measured pEW values for the Si II A4000 
feature within 5 d of maximum brightness. The pEWs are 
plotted against xi, c, and Hubble residuals corrected only 
for colour (for SNe la that are used when constructing the 
Hubble diagram). 

The pEW of Si II A4000 is highly correlated with xi 
(Figure |9j top plot) with a correlation coefficient of —0.87. 
This is in agreement with many previous studies and is ac- 
tually a stronger correlation than has been seen before (e.g., 
Arsenijevic et al.|[2008l [Walker et al.|[20TT| [Blondin et al.| 



2011[ [Nordin et ai.||2011b||Chotard et al.||2011| ). The leasT 
squares linear fit to the data is shown as the solid line in the 
top plot of Figure [9| (the dashed lines are the standard error 
of the fit). Note that the Ia-99aa objects in the figure also 
lie along the relationship. 

Arsenijevic et al.| ( [2008[ ) plotted xi versus the pEW of 
Si II A4000 and coded each low-redshift point based on its 
Benetti type; they found that the FAINT objects fell below 
the linear relationship (i.e., they had smaller than expected 
pEW values). In BSNIP II it was shown that FAINT (and 
similarly, Ia-91bg) objects have, if anything, larger than av- 
erage pEW values (especially for the Si II A4000 feature). 
Figure [9| shows no Ia-91bg objects since SALT/2 is unable 
to fit that spectral subtype. However, when the BSNIP val- 
ues of xi are plotted against the pEW of Si II A4000 and 
coded by Benetti type, the two FAINT objects fall at the 
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Figure 9. The pEW of the Si II A4000 feature versus SALT2 
light-curve width parameter xi (top), SALT2 colour parameter c 
(middle), and Hubble residuals corrected for colour only {bottom). 
Colours and shapes of data points are the same as in Figure [5| 
In the top plot, the solid line is the linear least-squares fit and 
the dashed lines are the standard error of the fit. In the middle 
plot, the median uncertainty in both directions is shown in the 
lower-right corner. 
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upper-left end of the linear correlation. The cause of this 
discrepancy between the two studies is unclear. 

The middle plot of Figure [9] shows no real evidence that 
the pEW of the Si II A4000 feature is correlated with c. The 
correlation coefficient we find for all of the objects is 0.13, 
and when removing objects with c > 0.5, the coefficient 
only increases to 0.15. This is slightly smaller than what was 
found by Blondin et al. j^2011J, an d significantly smaller than 
what was found by Nordin et al. (2011al. While the former 



claim no observed correlation, the latter do claim that the 
pEW of Si II A4000 is correlated with c. 

The bottom panel of Figure [9] shows the Hubble residu- 
als when corrected only for SALT2 colour versus the pEW of 
Si II A4000. Blondin et al. ( 2011| saw a relatively weak cor- 
relation between these parameters and found that a distance 
model involving c and the pEW of the Si II A4000 feature 
led to a "marginal improvement" over the standard {x\,c) 
model. We find a strong correlation (coefficient of 0.81), and 
the (c. Si II A4000 pEW) model performs nearly as well as 
the {xi,c) model (A^^^^c = 0.012 ± 0.036). T hus, the BSNIP 
data are in agreement with the finding of [Blondin et al.| 
(2011[) that the pEW of the Si II A4000 feature is essentially 



a replacement for the x\ parameter and is an accurate mea- 
surement of light-curve width. Figure [lO] shows the Hubble 
diagram residuals for the (c. Si II A4000 pEW) model and the 
colour-corrected-only model verus redshift, with the WRMS 
for each model as the grey band. 

Interestingly, we found that a model involving c, x\, 
and the pEW of the Si II A4000 feature actually leads to 
a ~ 10 per cent decrease in WRMS and a ~ 28 per cent 



decrease in cr 



prcd 



For this model, 



-0.026 ± 0.15, 



which implies that the improvement has a significance of 
about 1.8(7. The correlation between pEW of Si II A4000 and 
Hubble residuals corrected for colour and light-curve width 
is only 0.20, and thus perhaps the combination of c, xi, and 
pEW of the Si II A4000 feature is not actually adding much 
new information. However, the utility of a model including 
all three of these parameters should be investigated with 
other SN samples. Figure [TT] contains Hubble residuals for 
the {x\,c, Si II A4000 pEW) model as well as the standard 
(xi, c) model (using the same set of objects) versus redshift. 
Also shown, as the grey band, is the WRMS for each model. 

While the above investigation focused on the SALT2 
light-curve fitter, we can investigate correlations between the 
pEW of the Si II A4000 feature and photometric parameters 
from MLCS2k2. Figure shows the 63 BSNIP SNe which 
have MLCS2k2 fits and measured pEW values for the Si II 
A4000 feature within 5 d of maximum brightness. The pEWs 
are plotted against A and Av- 

The correlation coefficient between pEW of Si II A4000 
and A is 0.80, which is larger than previously observed 
(Nordin et al. 2011b I and implies that these parameters 



are highly correlated. This is expected based on the high 
degree of correlation between the pEW of Si II A4000 and 
x\ since both A and x\ are measurements of the width of 
SN la light curves. In Figures [12] and [9] the Ia-99aa objects 
lie at the extreme low-pEW end of the relationship, though 
perhaps they are slightly systematically below the trend in 
the top panel of Figure [12] On the other hand, the Ia-91bg 
objects in the top panel of Figure [12] fall significantly below 
the linear trend. 
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Figure 10. Hubble diagram residuals versus zcMB for the (c. Si II 
A4000 pEW) model {top) and the colour-corrected-only model 
{bottom). The grey band is the WRMS for each model. Colours 
and shapes of data points are the same as in Figure [5] 



relation between the pEW of Si II A4000 and Ay, even when 
objects with Ay > 0.5 mag are removed (correlation coeffi- 
cients of < 0.13 in both cases using the BSNIP data). In all 
of the plots presented in this section, HV and la-norm ob- 
jects overlap significantly. This is expected since in BSNIP II 
it was shown that the pEW of the Si II A4000 feature is ex- 
tremely similar for these two subclasses. We also note that 
the plot of Ami5(B) versus the pEW of Si II A4000 looks 
nearly identical to the top panel of Figure[T2]and has a larger 
correlation coefficient of 0.87. 



4.4.2 Fe II and Mg II 



Nordin et al. (2011b) found that the ApEW of Fe II within 



As in Nordin et al. (2011b I, there is no significant cor- 



3 d of maximum brightness is well correlated with SALT 
colour. In Figure [13] we present the 63 SNe in the BSNIP 
sample which have SALT2 fits and a measurement of the 
pEW of the Fe II complex within 5 d of maximum bright- 
ness. The two observables plotted have a correlation coeffi- 
cient of 0.49, increasing slightly for objects with c < 0.5. If, 
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Figure 11. Hubble diagram residuals versus zqmb for the (xi, c, 
Si II A4000 pEW) model (top) and the standard {xi,c) model 
{bottom). The grey band is the WRMS for each model. Colours 
and shapes of data points are the same as in Figure [s] 



however, only spectra within 3 d of maximum brightness are 
used, the correlation increases slightly but the sample size 
decreases by nearly one-quarter. The strength of this cor- 
relation is slightly higher than that found by |Nordin et al.] 
(2011b I, though we point out that they used SALT colour 



while we use SALT2 colour, and they used ApEW (i.e., the 
difference between the measured pEW and the average pEW 
evolution) while we use the actual measured pEW. Ia-99aa 
objects appear to have typical values of both the pEW of 
Fe II and c, though there are a very small number of objects 
of this spectral subtype. On the other hand, HV SNe seem to 
be both redder and have larger pEWs, but as before, there 
is significant overlap with la-norm objects as well. 

The pEW of the Mg II complex is also correlated with 
c, as seen in Figure [141 In that figure, 64 SNe within 5 d of 
maximum brightness are plotted, and they have a correla- 
tion coefficient of 0.56 (for objects with c < 0.5). Once again, 
when using only spectra within 3 d of maximum brightness 
the correlation becomes slightly stronger at the expense of 
significant decrease in the sample size. This correlation has 



Figure 12. The pEW of the Si II A4000 feature versus MLCS2k2 
light-curve width parameter A (top) and MLCS2k2 reddening 
parameter Ay (bottom). Colours and shapes of data points are 
the same as in Figure[5] In the top plot, the solid line is the linear 
least-squares fit and the dashed lines are the standard error of the 
fit. 



been observed previously, though at slightly lower signifi- 
cance and with almost one-third the number of low-z SNe 
( Walker et al.|20TT |. The data presented in Figure 14 match 
well with what was shown by Walker et al. (2011 1 for low-z 



objects, but their high-z sample spans a much larger range 
of Mg II pEW values and shows nearly no correlation with 
c. There seems to be effectively no difference between the 
various spectroscopic subtypes in this parameter space. 

The significant correlations between c and both the 
pEWs of Mg II and Fe II are promising. Measuring a 
pEW of a broad feature in a single spectrum near maxi- 
mum brightness is much simpler than obtaining photomet- 
ric data for a full light curve that needs to then be mod- 
eled by a light-curve fitter (such as SALT2). Therefore, as 
with all other spectral measurements discussed in this work, 
we constructed Hubble diagrams using the pEW of Mg II 
and Fe II and models of the form shown in Equations |4]- 
[7] All but one of these models performed worse than the 
standard (a;i,c) model. The model including xi, c, and the 
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Figure 13. The pEW of the Fe II complex versus the SALT2 
colour parameter c. Colours and shapes of data points are the 
same as in Figure [s] The median uncertainty in both directions 
is shown in the lower-right corner. 




Figure 14. The pEW of the Mg II complex versus the SALT2 
colour parameter c. Colours and shapes of data points are the 
same as in Figure [5] The median uncertainty in both directions 
is shown in the lower-right corner. 



pEW of Mg II was only as accurate as the standard model 
(A^i,c = -0.004 ±0.004). However, using the pEW of Mg II 
or Fe II as a replacement for c or in addition to c (along with 
xi) is a tantalizing possibility that should be explored using 
future datasets. 

Furthermore, we attempt two Hubble diagrams using 
no light-curve information whatsoever. One uses only the 
pEWs of the Si II A4000 feature and the Mg II complex, and 
the other uses only the pEWs of the Si II A4000 feature and 
the Fe II complex. The idea is that the pEW of Si II A4000 
is a good proxy for xi , and the pEWs of the Mg II and Fe II 
features are reasonably good proxies for c. These Hubble 
diagrams included only a subset of the data used in the 
flux-ratio study (see Section [4T0| since they needed to have 
well-measured pEWs. While both the Mg II and Fe II models 



had quite low WMRS values (0.274 and 0.297, respectively), 
they were not as low as the WRMS values using the standard 
(xi,c) model (0.118 and 0.121, respectively). 

44.3 SII "W" 

As discussed in Section |4.3[ the depth of the bluer absorp- 
tion of the S II "W" feature relative to the pseudo-continuum 



was shown by Blondin et al. (20111 to decrease the scatter 



of Hubble residuals by about 10 per cent. In that section 
we showed that the relative depths of the redder absorption 
of the S II "W" in the BSNIP data were marginally corre- 
lated with xi, opposite to what was seen by [Blondin et alT] 
(12011 1. However, both studies agree that c and the colour- 
and width-corrected Hubble residual are uncorrelated with 
the relative depth of the S II "W." 

As discussed in BSNIP II, the relative depth of a spec- 
tral feature relies on a spline fit to the spectra and can fairly 
easily be contaminated by local noise. The pEW, however, 
is less prone to this type of contamination, relies only on the 
definition of the pseudo-continuum (and not any additional 
fit to the data), and often contains the same information 
as the relative depth. For these reasons the pEW values 
were used in favor of the a values in the analysis performed 
in BSNIP II. Furthermore, both I Blondin et H] ( |20TT| ) and 
BSNIP II measure the pEW of the entire S II "W" feature, 
and thus pEW values are a more fair comparison between 
the two studies than are a values. 

We find that c is uncorrelated with the pEW of the S II 
"W" (correlation coefficient of 0.15 for objects with c < 0.5). 
In contrast with the relative depth of the S II "W," however, 
xi is also uncorrelated with the pEW of the S II "W" (corre- 
lation coefficient of —0.16). This is significantly weaker than 
the correlation between xi and the relative depth of this fea- 
ture found in Section [4.3[ Equations |4||7] were used to create 
Hubble diagrams involving the pEW of the S II "W," but 
none of these models led to an improvement in the WRMS. 

4.4.4 Si II X5972 

In BSNIP II it was shown that the pEW of the Si II A5972 
feature correlated well with the spectral luminosity indica- 
tor 5R(Si II) (see Section |4.5| for more information on this 
parameter). Thus, one might expect this pEW to be an ac- 
curate luminosity indicator as well, and in fact evidence for 
a correlation between the pEW of Si II A5972 and both xi 



and Ami5(_B) has been seen in previous work (Nordin et al. 
2011b[|Hachinger et al.|2006j) 



Figure |15| shows the 55 SNe which have a SALT2 fit 
as well as a measured pEW for the Si II A5972 feature. We 
find a correlation coefficient of —0.67, which is stronger than 
what was found by |Nordin et al.| ( |2011b[ ). As in the relation- 
ship between xi and the pEW of Si II A4000, the Ia-99aa 
objects appear to follow the relation. 

Similarly, a strong linear correlation has been observed 
betwe en the pEW of Si II A5972 and Ami5(_B) (Hachinger] 



et al.j|2006). This relationship is also found in the BSNIP 



data, as shown in Figure [16) there are 62 SNe, the param- 
eters have a correlation coefficient of 0.67, and as with the 
relationship with xi, Ia-99aa objects occupy the bottom of 
the correlation while the la-norm and HV objects are highly 
overlapping. 
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Figure 15. The pEW of Si II A5972 versus the SALT2 hght-curve 
width parameter xi. Colours and shapes of data points are the 
same as in Figure |5] 
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Figure 16. The pEW of the Si II A5972 feature versus Ami5(B). 
Colours and shapes of data points are the same as in Figure [5] 



Hachinger et al. ( 2006 1 denote the Benetti type of each 
object on their plot of the pEW of Si II A5972 versus 
Ami5(B) and note that FAINT objects are found at the 
top of the correlation while LVG SNe are found at the bot- 
tom (with most HVG objects occupying the middle of the 
trend). The BSNIP data show a similar behaviour for the 
FAINT objects (again, much like the Ia-91bg objects in Fig- 
ure 161; however, our data show no differentiation between 



the LVG and HVG objects in this parameter space. 

Finally, we note that when plotting the pEW of the Si II 
A5972 feature against the MLCS2k2 A parameter, nearly the 
exact same results are seen as those in Figure [16] As with 
the S II "W," no distance model utilising the pEW of Si II 
A5972 led to an improvement in the Hubble residuals. 



4.4.5 St 11X6355 

Much like the pEW of Si II A5972, the pEW of Si II A6355 
has been seen to correlate marginally well with xi and 
to separate various spectral subtypes when compared to 
Ami5(B) ( [Nordin et "aL]|2011b| [Hachinger et al.||2006"l re- 
spectively) . 

In Figure [l7[ the 66 SNe with SALT2 fits and pEW 
values for the Si II A6355 feature are shown. The pEW values 
are plotted against xi, c, and Hubble residuals corrected for 
light-curve width and colour (for objects that are part of 
the Hubble diagram). A correlation coefficient of —0.59 is 
calculated for the top panel. This is consistent with what 



was observed in the data studied by Nordin et al. (2011b I 



As in the previous two relationships between xi and the 
pEW of Si II features, the Ia-99aa objects follow the linear 
relation. 

We find that c is somewhat correlated with the pEW 
of Si II A6355 (correlation coefficient of 0.41) for objects 
with c < 0.5. This pEW is even less correlated with Hubble 
residuals corrected for xi and c (correlation coefficient of 
0.23). A distance model which includes xi, c, and the pEW 
of the Si II A6355 feature leads to a 4 per cent decrease in 
WRMS, a 6 per cent decrease in o-pred, and is significant at 
the 1.2(7 level. So while this is technically an improvement 
over the standard {x^ , c) model, it may not actually be very 
helpful. 

Plotting the pEW of Si II A6355 versus Ami5(B), 
Hachinger et al.[ ( [2006[ ) are able to separate FAINT, LVG, 
and HVG objects relatively accurately. This is also seen, 
though at a lower significance, in the BNSIP data. In Fig- 
ure [T8[ we plot 80 SNe; the parameters have a correlation 
coefficient of 0.52, but the Ia-91bg, Ia-99aa, and the lone la- 
91T objects are all reasonably well separated from the bulk 
of the SNe. There is even some evidence for a difference be- 
tween HV and la-norm objects in this parameter space. As 



mentioned above, Hachinger et al. 1 2006 1 denote the Benetti 



type of each object on their plot of the pEW of Si II A6355 
versus Amis{B) and state that the three subtypes are well 
separated. Again the BSNIP data support this conclusion, 
but at a weaker significance. FAINT objects are found in the 
same part of parameter space as the Ia-91bg objects, while 
HVG and HV SNe tend to occupy a different part of param- 
eter space compared with the LVG and Ia-norm/91T/99aa 
objects. However, there is quite a lot of overlap among all 
of the non-FAINT (and non-Ia-91bg) SNe. 

Figure |19[ presents the pEW of the Si II A6355 feature 
plotted against (B — V)^^^. The top panel shows all 78 SNe 
from the BSNIP dataset for which both of these values have 
been measured and the bottom panel displays a close-up 
of objects with (B — V)^^^ < 0.319 mag. The linear least- 
squares fit to all of the points is represented by the solid line, 
and the fit to SNe with {B — V^)^^;^,^ < 0.319 mag is shown 
by the dotted line. 

The correlations for the full sample and the less red- 
dened sample are weak (coefficients 0.12 and 0.16, respec- 
tively). Qualitatively, this matches the work of [Foley et al.[ 
poTi| ), even though they plot intrinsic colours and we plot 
observed colours. The pEW of the Si II A6355 feature is less 
correlated with observed {B — V)^^^ than the velocity near 
maximum brightness of that same spectral feature. 
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Figure 17. The pEW of the Si II A6355 feature versus SALT2 
light-curve width parameter xi {top), SALT2 colour parameter 
c {middle), and Hubble residuals corrected for light-curve width 
and colour {bottom). Colours and shapes of data points are the 
same as in Figure |5] 



Figure 18. The pEW of Si II A6355 versus Ami5(B). Colours 
and shapes of data points are the same as in Figure [5] 



44.6 Ca II and O I 



Much like [Hachinger et aL\ ( |2006[ ), we searched for possi- 
ble correlations between the pEW of each spectral feature 
investigated and various photometric parameters. Many of 
the strongest and most interesting of these possible corre- 
lations have been discussed in the preceding sections. For 
Ca II H&K as well as the O I triplet, no pairs of pEWs 
and photometric parameters were found to have correlation 
coefficients > 0.4. 

However, the pEW of the Ca II near-IR triplet is found 
to correlate with A (with a correlation coefficient of 0.53) 
and with c (also a correlation coefficient of 0.53 for SNe 
with c < 0.5). Both correlations are marginal at best; how- 
ever, this spectral region has been studied very little in the 
past. One of the strengths of the BSNIP data are that the 
average wavelength coverage (3300-10,400 A, BSNIP I) is 
significantly wider than that of most other SN la spectral 
datasets. For example, one of the largest previously pub- 
lished SN la spectral dat asets had an average wavelength 
coverage of 3700-7400 A ( [Matheson et al.||2008[ ). Thus, the 
Ca II H&K feature, the O I triplet, and the Ca II near-IR 
triplet have been ignored almost entirely in past spectral 
analyses like the one presented here. 

We once again constructed Hubble diagrams from 
Equations |4}|7] and using the pEW values of Ca II H&K, 
the O I triplet, and the Ca II near-IR triplet. All but 
two models were significantly worse at measuring distances 
than the standard {xi,c) model. The {xi,c, Ca II H&K 
pEW) and (a;i,c, O I triplet pEW) models both slightly 
decreased the WRMS, but at almost imperceptible levels 
(A^i,^ = -0.0023 ± 0.0096 and A^i,^ = -0.0063 ± 0.0148, 
respectively) . 



4.5 The Si II Ratio 

Historically, one of the first spectral luminosity indicators 



investigated was the Si II ratio, 5R(Si II), defined by Nugent 
|et al.| ( |1995[ ) as the ratio of the depth of the Si II A5972 
feature to the depth of the Si II A6355 feature. IHachinger] 
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Figure 19. The pEW of the Si II A6355 feature ver- 
sus {B — V)^^-^ (top) and a close-up view of objects with 
{B — V)^^^ < 0.319 mag {bottom). Colours and shapes of data 
points are the same as in Figure [5] The solid line is the fit to all 
of the data while the dotted line is the fit only to objects with 
(B-V) max 

< 0.319 mag. 



et al.| (|2006| redefined the Si II ratio as the pEW of Si II 
A5972 divided by the pEW of Si II A6355. In BSNIP II it 
was shown that these are nearly equivalent definitions, so 
in order to be consistent with that work we define the Si II 
ratio for the present study to be 



K (Si II) 



pEW (Si II A5972) 
pEW (Si II A6355) ' 



(12) 



The Si II ratio has been shown to correlate with max- 
imum absolute B-band magnitude and Ami5(B), which is 
why it has been used as a spectral luminosity indicator (e.g., 
Nugent et aLl [19951 [Benetti et "aLj [20051 [Hachinger et al.| 
2006[). Figure (20[ shows 62 SNe la with both Ami5(B) and 



5R(Si II). The data are correlated with a correlation coeffi- 
cient of 0.47. Ia-91bg objects appear at the upper right of 
the plot and form a continuous relationship with the la-norm 
objects. Ia-99aa objects appear to lie above the main trend 
(though there are only a handful of these SNe in the figure), 
while HV objects lie below the main trend. When removing 



Figure 20. The Si II ratio versus Ami5{B). Colours and shapes 
of data points are the same as in Figure [sl 



the Ia-99aa objects, the correlation increases slightly (coef- 
ficient of 0.54). 

If we instead tag each data point in Figure [20l by its 
Benetti type, we find that the FAINT objects are found in 
the upper right of the main trend (similar to the Ia-91bg ob- 
jects and as seen in previous studies, e.g., [Benetti et al.|2005| 
Hachinger et al. 2006 1. The LVG and HVG objects are found 



in the lower-left portion of the plot with a significant level 
of overlap between the two subclasses. This differs from pre- 
vious work, where there have been claims that LVG objects 
have larger Si II ratios and lie above the main trend ( [Benettil 
et al.|2005 Hachinger et al.|2006[ ), though these studies and 
the BSNIP data both observe larger scatter in 5R(Si II) val- 
ues in the lower Ami5(B) objects. When removing the LVG 
objects, the correlation is effectively unchanged. Comparing 
5R(Si II) to the MLCS2k2 A parameter yields similar results 
to those seen in Figure [20[ 

The BSNIP distribution of Ami5{B) values, while not 
evenly distributed, is more continuous than in previous 
studies similar to the present one. For example, the data 
in Hachinger et al. (20061 contained only one object with 



Ami5(_B) between 1.5 and 1.7 mag, while the BSNIP data 
have 7 objects in that range. A more continuous distribu- 
tion of Ami5(i3) values, combined with the spectroscopic 
subclasses presented in Figure [20[ complicates the relatively 
simplistic view that underpins the basic Phillips relation. For 
0.95 < Amis{B) < 1.0 mag there are Ia-99aa, la-norm, and 
HV objects. On the other end of the Amis{B) distribution, 
between ~1.75 and ~1.95 mag there are Ia-91bg, la-norm, 
and again HV objects. Thus, for a given light-curve width 
(or decline rate) , there exist SNe la of significantly different 
subclasses. 

As discussed in BSNIP II, objects tagged by SNID as 
Ia-91bg or Ia-99aa are the most spectroscopically peculiar 
objects and probably only represent the extreme ends of a 
continuous distribution of spectra. If true, this means that 
the most spectroscopically peculiar objects may not have 
the most extreme light curves. The reverse may also be true, 
namely that the SNe with the most extreme light curves may 
not be the most spectroscopically peculiar. This is further 
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supported by the relatively wide scatter in the main trend 



of Figure 20 At any value of Ami5(_B) (with a significant 



number of objects) there is a broad range in 5R(Si II) values. 

In BSNIP II, it was pointed out that the relative 
strength of the two Si II features that go into calculating 
5R(Si II) is fairly robust at differentiating between the various 
"SNID types" and "Wang types." Thus, from a spectrum, 
either using SNID or the pEWs of Si II features, one may 
declare an object to be Ia-91bg or Ia-99aa, whereas based 
on the light curve of the same object it might be considered 
relatively normal. The significant amount of scatter in the 
correlation between the Si II ratio and Ami5(_B) also cau- 
tions one against simply measuring 5R(Si II) from a single 
spectrum and then using that value and a fit to the data in 



Figure 20 to calculate a Ami5(_B) value. 



In Figure pi] we show the 51 BSNIP SNe which have 
SALT2 fits and Si II ratios within 5 d of maximum bright- 
ness. 5R(Si II) is plotted against x\, c, and Hubble residuals 
corrected for colour only (for SNe which are used in the 
Hubble diagram). 

From the BSNIP data we find that the Si II ratio is only 
marginally correlated with xi (top plot of Figure |21[ ) with 
a correlation coefficient of —0.42. This is a weaker relation- 
ship than what has been found before ( Blondin et al.|20"TT I . 
The Ia-99aa objects appear to be above the main trend and 
most of the HV objects seem to be below it. The middle 
plot of Figure |2l] shows no evidence for a correlation be- 
tween 5R(Si II) and c, with a correlation coefficient of 0.07 
for objects with c < 0.5. 

The bottom plot of Figure [2l] shows a possible correla- 
tion between the Si II ratio and Hubble residuals corrected 
for colour only (coefficient of 0.34), and it is again signif- 
icantly weaker than what has been found before (Blondin 
et al.|20TT l. In fact, Blondin et al. (2011 1 go so far as to say 
that 5R(Si II) acts as a replacement for x\, but the BSNIP 
data do not support such a claim. Using the current sample, 
the best model which includes the Si II ratio also includes 
both x\_ and c, but it is only about as accurate as the stan- 
dard (a;i,c) model (A^,,e = 0.0199 ± 0.0199). 



Interestingly, Blondin et al. (20111 found that the sub- 
luminous (but la-norm) SN 2000k is a 2a outlier in their 
plot of K(Si II) versus c, but part of the main correlation of 
5R(Si II) versus x\. This object is in the BSNIP dataset and, 
while we agree that it is subluminous and spectroscopically 
normal, it is not a significant outlier in any of the three plots 
in Figure [21] 



4.6 The Ca II Ratio 



The Ca II ratio was defined by [Nugent et al.| ( [1995[ ) as the 
ratio of the flux at the red edge of the Ca II H&K feature 
to the flux at the blue edge of that feature. In the notation 
from BSNIP II this is 



5R (Ca II) 



Fr (Ca II HfcK) 
Fb (Ca II H&K) ■ 



(13) 



Like the Si II ratio, it has been found to correlate with max- 



imum absolute _B-band magnitude ( Nugent et al.[[1995 K 

In BSNIP II it was shown that the Ca II ratio and the 
Si II ratio are uncorrelated, even though both of them have 
been used as spectral luminosity indicators. In Section [4.5| 
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Figure 21. The Si II ratio versus SALT2 light-curve width pa- 
rameter xi (top), SALT2 colour parameter c {middle), and Hubble 
residuals corrected for colour only {bottom). Colours and shapes 
of data points are the same as in Figure[5] In the middle plot, the 
median uncertainty in both directions is shown in the upper-right 
corner. 
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Figure 22. The Ca II ratio versus Ami5(_B). Colours and shapes 
of data points are the same as in Figure [5] The solid line is the 
linear least-squares fit and the dashed lines are the standard error 
of the fit. 



it was shown that 5R(Si II) is correlated with Ami5{B). Fig- 
illustrates that 5R(Ca II) is correlated with Ami5(_B) 



22 



as well (correlation coefficient of 0.67). The plot contains 65 
SNe; the least-squares linear fit to the data is shown as the 
solid line in Figure [22] and the dashed lines are the standard 
error. 

Interestingly, the HV objects seem to occupy a rela- 
tively narrow region of parameter space that is surrounded 
on all sides by mainly la-norm SNe. Ia-99aa objects mostly 
make up the lowest end of the linear trend, while one of the 
two Ia-91bg objects in the plot perhaps does not follow the 
main relationship. Unsurprisingly, comparing A to 5R(Ca II) 
results in the same trends seen in Figure [22] 

Figure |23] displays the 64 BSNIP SNe which have 
SALT2 fits as well as Ca II ratios within 5 d of maximum 
brightness. 5R(Ca II) is plotted against xi, c, and Hubble 
residuals corrected for light-curve width and colour (for ob- 
jects used when constructing the Hubble diagram). 

The Ca II ratio appears to be well anticorrelated with 
xi (correlation coefficient of —0.58) and correlated with c 
(coefficient of 0.48). The bottom plot of Figure 



23 shows that 



the Xi and c corrected Hubble residuals and K(Ca II) are not 
well correlated with each other (coefficient of 0.24). However, 
a model that uses xi, c, and the Ca II ratio decreases the 
WRMS by ~6 per cent and the cTpred by ~33 per cent. The 
significance of this improvement is only at the l.lcr level 
(A^i,c = -0.0207 ±0.0191), however. 
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4.7 The "SiS" Ratio 

Analogous to the Ca II ratio, the "SiS ratio" was introduced 



by Bongard et al. ( 2006 1 as the ratio of the flux at the red 



edge of the S II "W" feature to the flux at the red edge of 
the Si II A6355 feature. In the notation used in BSNIP II 
this is 



K(SiS) 



Fr (S II "W") 
Fr (Si II A6355) ' 



(14) 



Figure 23. The Ca II ratio versus SALT2 light-curve width pa- 
rameter xi (top), SALT2 colour parameter c {middle), and Hub- 
ble residuals corrected for light-curve width and colour {bottom). 
Colours and shapes of data points are the same as in Figure [5] 
In the middle plot, the median uncertainty in both directions is 
shown in the lower-right corner. 
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Figure 24. The SiS ratio versus Ami5(iJ). Colours and shapes of 
data points are the same as in Figure[5] The correlation coefficient 
is -0.54. 



Figure 25. The SiS ratio versus c. Colours and shapes of data 
points are the same as in Figure [5] The median uncertainty in 
both directions is shown in the upper-right corner. 



In a sample of 8 SNe, K(SiS) has been seen to correlate with 
maximum absolute _B-band magnitude in the same way as 
K(Call) ( [Bongard et al.|2006[ ). 

The SiS ratio and the Si II ratio were found to be only 



marginally correlated in BSNIP II, and Figure 24 (which 



contains 72 SNe) shows that SR(SiS) is anticorrelated with 
Ami5(B) (correlation coefficient of —0.54). Note that this 
relationship is in the opposite sense of the one between 
the Ca II ratio and /S.mi^{B). Here Ia-99aa/91T objects lie 
above the main relationship while the la-norm and HV SNe 
are well mixed. Once again, comparing A or x\ to 5R(SiS) 
yields similar results to what is seen in Figure |24] 

The SiS ratio appears to be well correlated with c, when 
including the most reddened objects. In Figure [25] are 71 
SNe, and the correlation coefficient is —0.50. However, if one 
removes the most highly reddened objects with c > 0.5, the 
correlation weakens slightly to —0.47. No distance model in- 
volving the SiS ratio is more accurate than the {xi , c) model. 
However, when 5R(SiS) is combined with just c or both xi 
and c, the accuracy is on par with the standard (xi, c) model 
(A^i,c = 0.0209 ± 0.0223 and A^, = 0.0074 ± 0.0105, re- 



spectively). We also find (in Section 4.10 1 that out of 17,822 
flux ratios combined with c, the most accurate distances are 
calculated using flux ratios that are effectively the SiS ratio, 
and that it is nearly as accurate as using the standard {xi, c) 
model. Finally, we note that |Blondin et al.| ( |2011 | found that 
5R(SiS) performs significantly worse than the usual (xi,c) 
model. 



4.8 The "SSi" Ratio 

Yet another possible spectroscopic luminosity indicator is 
the ratio of th e pEW of the S II "W " to that of the Si II 
A5972 feature ( [Hachinger et al.||2006[ ). This SSi ratio is de- 
fined in BSNIP II as 



K(S,Si) 



pEW(S II "W") 
pEW (Si II A5972) ' 



CO 




1.5 

Am, 5(B) (mag) 

Figure 26. The SSi ratio versus Ami5(B). Colours and shapes 
of data points are the same as in Figure [s] 



Hachinger et al. (20061 found that the SSi ratio is linearly 



anticorrelated with Amis (which is opposite to the relation- 
ship between K(Si II) and Amis). The analysis in BSNIP II 
seemed to confirm this observation by showing that the SSi 
ratio was strongly anticorrelated (nonlinearly) with the Si II 
ratio. 

5ft(S,Si) is plotted agai nst Amis (5) for 59 SNe in Fig- 
ure 26 The results of Hachinger et al. ( 2006 1 and the specu- 



(15) 



lation in BSNIP II are confirmed: the SSi ratio is strongly an- 
ticorrelated with Ami5(_B) (correlation coefficient of —0.65). 
Here, the Ia-91bg and Ia-99aa objects follow the main trend 
and are found at the lower and upper ends of the correla- 
tion, respectively. There are only a few HV objects in Fig- 
ure |26[ but there is some evidence that they have larger 
than average 5R(S,Si) values (which was also seen in Fig. 16 
of BSNIP II). 

Plots of A and xi versus 5R(S,Si) display trends like that 
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Figure 27. The SiFe ratio versus Ami5{B). Colours and shapes 
of data points are the same as in Figure [5] The solid line is the 
linear least-squares fit and the dashed lines are the standard error 
of the fit. 



of Figure [26] However, the Ia-99aa objects fall off of the main 
correlation in both of these parameter spaces. In both cases 
these SNe have lower 5R(S,Si) values than one would expect 
from the main correlation. 



4.9 The "SiFe" Ratio 

Analogous to the SSi ratio, the "SiFe ratio" was defined 
as the ratio of the pEW of the Si II A5972 feature to that 
of the Fe II complex, and it was shown to be an accurate 
spectroscopic luminosity indicator ( jHachinger et al^|2006[ ). 
In BSNIP II, K(Si,Fe) was defined as 



K (Si,Fe) 



pEW (Si II A5972) 
pEW (Fe II) ' 



(16) 



and found to be strongly correlated with the Si II ratio. 

We plot SR(Si,Fe) versus AmisjB ) for 5 3 SNe in Fig- 
ure |27] The results of [Hachinger et ah] ( |2006[ ) and the spec- 
ulation in BSNIP II are again confirmed: the SiFe ratio is 
strongly (linearly) correlated with Ami5(i3), with a corre- 
lation coefficient of 0.60. The solid line in the figure is the 
linear least-squares fit and the dashed lines are the standard 
error of the fit. 

Similar to 5R(Si II), Ia-99aa objects are found at the 
lowest end of the linear trend while the Ia-91bg objects in 
the plot appear to be above the main relationship. In Fig- 
ure [27| there are only a few HV SNe, but they appear to 
have smaller than average 5R(Si,Fe) values (which was also 
seen in Fig. 17 of BSNIP II). When comparing A and xi to 
5R(Si,Fe), the basic trend seen in Figure [27] is recovered, but 
with larger scatter (even though the correlation coefficients 
are similar). 

4.10 Arbitrary Flux Ratios 



Bailey et al. ( 2009 \ found that by using ratios of fluxes from 



as TZ(\y/\x) = F{\y) / F{\x), where Xy and A^, are the rest- 
frame central wavelengths of given bins[^ The spectra are 
forced to cover a wavelength range of exactly 3500-8500 A 
and are binned into 134 equal-sized (in In A) bins (corre- 
sponding to 2000 km s~^ per bin). The data are also dered- 
shifted and dereddened using the redshift and reddening val- 
ues presented in Table 1 of BSNIP I and assuming that the 



extinction follows the Cardelli et al. ( 1989|) extincti on law 



modified by O'Donnell (19941. As in Bailey et al. (20091 



and Blondin et al. (2011 1, a colour-corrected version of this 



fiux ratio {TV^ {Xy / \x)) is also used, and it is defined as the 
ratio of fluxes as measured from a spectrum that has been 
corrected for SALT2 c using the colour law from |Guy et al.| 
( pOOT] ). We use these colour-corrected flux ratios when test- 
ing models that also adopt the SALT2 colour parameter (i.e.. 
Equations [6] and [t]) . 

As with the rest of the current study, we only investi- 
gate spectra within 5 d of maximum brightness since it was 
shown in BSNIP II that the spectra do not evolve signifi- 
cantly during these epochs. Also, as mentioned above, we do 
not use only spectra within 2.5 d of maximum (as has been 
done previously. Bailey et al.|2009 Blondin et al.|2(5TT l be- 
cause the significance of our results would be weakened due 
to the smaller number of objects. Since the average spec- 
trum in BSNIP extends to 3300 A (BSNIP I), we perform 
the current fiux-ratio analysis with the requirement that all 
spectra cover a wavelength range of 3300-8500 A. No ratios 
involving wavelengths below 3500 A are found to decrease 
the WRMS significantly. We also vary the binning of the 
spectra used in the fiux-ratio analysis and investigate data 
with bin sizes of 4000, 8000, and 10,000 km s~\ We find 
that as the spectra are binned more, the WRMS values in- 
crease for all ratios. This can be explained by the idea that 
larger bins will "blend" wavelength bins of fiux ratios that 
decrease the Hubble residuals with ones that do not and will 
thus add "noise" into flux ratios. 

Since we are utilizing SALT2 fits and Hubble diagrams, 
we again require th at SNe have ZheUo > 0.01, c < 0.50, 
and reduced < 2. Blondin et al. (20111 also require that 



the absolute difference between B — V colour at maximum 
brightness derived from the spectrum and derived from the 
photometry be less than 0.1 mag. This is used as a proxy 
for their relative spectrophotometric accuracy. In BSNIP I 
it was shown (in Table 3) that the relative spectrophoto- 
metric accuracy is often < 0.1 mag for the BSNIP data. In 
fact, B — V colour is only inaccurate at the 0.1 mag level 
for the oldest [t > 20 d) and noisiest (S/N < 20) BSNIP 
spectra. Therefore, the spectra investigated here should all 
be spectrophotometrically accurate enough for the flux-ratio 
analysis. 

Of the data studied here, 62 objects have flux ratios 
calculated for the entire wavelength range mentioned above 
and reliable SALT2 flts that pass our Hubble diagram crite- 



ria (see Section 3.3 1. We randomly divide our sample into 9 



a single, binned SN la spectrum they could decrease the 
scatter in their Hubble diagrams. These ratios are defined 



^ This is similar to the definition of [Bailey et al] |2009[|, bu t 
is the reciprocal of the definition used by [Blondin et al.| ( |2011[ |. 
However, this only really matters for the plots of Aj, versus Xx in 
Figure p9| Thus, each panel in Figure 7 of [Blondin et al.[ | |20lT| l 
is the transpose of the panels in Figure |29| When using either 
definition, note that the first wavelength listed for a given TZ. is 
the numerator in the actual ratio of fluxes. 



BSNIP III: Spectral and Photometric Analysis 23 



=1 



0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 













— 1 











0.02 



0.04 0.06 

^CMB 



0.08 



0.10 



Figure 28. Hubble diagram residuals versus zcMB for the stan- 
dard (xi, c) model (Equation [sjl. The grey band is the WRMS for 
the model. Colours and shapes of data points are the same as in 
Figure [5] 

groups of 7 SNe and 1 group with 6 SNe when doing 10-fold 
CV. 



4.10.1 Flux-Ratio Results 

The "best" flux ratios for each model are chosen to be the 
ones with the lowest WRMS values. Ranking by other pa- 
rameters, such as the intrinsic prediction error (as used in 
[Blondin et al.|[2011[ ), yields different values for the best- 
performing flux ratios. However, since our main goal is to 
minimise the scatter in the Hubble diagram, and since the 
WRMS has a relatively straightforward interpretation, we 
rank the best flux ratios by their WRMS values. As dis- 
cussed in Section [3.3| for each model involving a flux ratio 
(Equations |4||7| we calculate, in addition to the WRMS, 
the intrinsic prediction error (o-prcd), the intrinsic correla- 
tion {pxi,c) of the residuals with residuals using the (a;i,c) 
model (Equation |8|, and the difference (A^^.c) in intrinsic 
prediction error with respect to the {xi , c) model and its sig- 
nificance. These parameters, along with the wavelengths, of 
the top 10 ratios for each model which includes a flux ratio 
(Equations |4||7| , are shown in Table [l] 

Also displayed in Table [l] is the WRMS and Uprcd of 
our benchmark {xi,c) model. Figure [28| shows the Hub- 
ble diagram residuals for this model versus redshift for the 
62 SNe la mentioned above. The grey band indicates the 
WRMS for the model. 



Figure 29 shows the WRMS (left column) and abso- 
lute Pearson correlation coefficient of the correction term 
(either 'yTZ or 'yTZ'^) with the uncorrected Hubble residuals 
(right column) for aU 17,822 (= 134 x 133) flux ratios in aU 
four models involving a flux ratio: TZ, {xi,TZ), {c,TV^), and 
{xi,c,TZ'^) (top to bottom, respectively). 

The left column of Figure [29] is a proxy for the overall 
scatter in the model while the right column indicates how 
much "new" information is gained by adding in the correc- 
tion term {'yTZ or 'yTZ'^). In general, a model may have a low 
WRMS (or (Tprod), meaning that the model is fltting the data 



well, but since the data have uncertainties associated with 
them, the model might be overfitting the data and actually 
end up fitting noise. 

One way to discern whether this is the case is to see 
how well the correction terms correlate with the uncorrected 
Hubble residuals. As described in Section [4.2.2| if the terms 
that do not contain Xc or c (i.e., 'yTZ or "/TZ'^) are well corre- 
lated with the uncorrected Hubble residuals, then the mea- 
sured observable is fitting information that actually exists in 
the data (as opposed to noise). However, a large correlation 
does not necessarily imply a good model. This is obvious (for 
example) in the top row of Figure 29 where most flux ra- 



tios have large WRMS values, including ones that also have 
quite high correlations between the correction term and the 
uncorrected Hubble residuals. 



4.10.2 Model 1: TZ Only 

Using only a flux ratio (Equation|4| leads to no improvement 
over the usual c) model (Equation[8|. In fact, this model 
seems to perform signiflcantly worse, as can be seen by the 
relatively large A^^^c values. The WRMS and ffpred of the 
"best" ratios are quite a bit larger than those of the {xi , c) 
model. This differs from the conclusion of previous work, 
which found that models using a flux ratio alone could per- 
form as well as the {xi , c) model f Baile y et al.|200"9 Blondin 



|et al.|20TT| . The best-performing ratio in the 7?.-only model, 
7?. (7770/3750), is not correlated with xi, but is highly cor- 
related with c (correlation coefficient of 0.81). This implies 
that 7?. (7770/3750) is effectively a colour indicator. 

The best ratio using this model seen by [Blondin et alT] 
( 2011[ ), TZ (6630/4400), was found in their study to be simi- 
larly correlated with SALT2 colour and to improve the Hub- 
ble diagram residuals over using the (xi, c) model (albeit at a 
low significance). We tested all of our flux-ratio models using 
a randomly selected subset of 26 SNe from the BSNIP sam- 
ple in order to match the number of objects used in |BlondTnl 
|et al.| ( |20lf| ). These mod els all yielded WRMS values simi- 
lar to what was found by Blondin et al. (2011 1, though the 



best-performing flux ratios were significantly different and 
the significance of the improvement was decreased (mostly 
due to the fact that the (xi, c) model performed much worse 
when using only 26 objects). Much like our analysis which 
only used spectra within 2.5 d of maximum brightness, the 
decrease in sample size was chiefiy responsible for the weak- 
ened significance of Hubble residual improvement. 

[Bailey et al.| ( |2009| l also saw an overall improvement 
when using their top ratio for this model, 7?. (6420/4430). 
The WRMS for both of these ratios using the BSNIP data is 
~0.35 mag, which is much larger than the WRMS values seen 
for our best ratios. When the BSNIP data are flt using K = 2 
(which is what is effectively done by Bailey et al. 20091, 
the WRMS of the best-performing ratio does decrease as 
compared to our 7^ = 10 run. However, the best-performing 
ratios involve wavelengths that are randomly scattered in 
wavelength space. Thus, it appears that this oversimplifled 
approach of splitting the data evenly into one "training" set 
and one "validation" set is decreasing the Hubble residuals 
by fitting to noise, and so these results are less reliable than 
the more statistically rigorous run using K = 10. 
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Figure 29. A map of the WRMS {left column) and the absolute Pearson correlation coefBcient of the correction term with the uncorrected 
Hubble residuals {right column) for all 17,822 flux ratios used in all four models involving a flux ratio (7?. only, (xi,7?.), (c, T?.'^), and 
(o-'i, c, T?-'^), top to bottom, respectively). 
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Table 1. Top 10 Flux Ratios for Each Model 



Rank 


Ay 




WRMS (mag) 


o-prcd (mag) 


Pxi ,c 






n 


1 


7770 


3750 





218 ±0 


027 


0.179 ±0 


023 


-0.11 


0.095 ±0.030 


(3.2(t) 


2 


7670 


3750 





223 ±0 


028 


0.189 ±0 


023 


0.09 


0.104 ± 0.030 


(3.5cr) 


3 


7720 


3750 





223 ±0 


027 


0.188 ±0 


023 


-0.05 


0.103 ±0.030 


(3.4cr) 


4 


7930 


3750 





226 ±0 


030 


0.187 ± 


023 


-0.10 


0.104 ± 0.031 


(3.4cr) 


5 


6990 


3750 





227 ± 


030 


0.193 ±0 


023 


-0.17 


0.111 ± 0.031 


(3.60-) 


6 


7880 


3750 





227 ± 


027 


0.190 ±0 


023 


-0.06 


0.106 ±0.030 


(3.5(t) 


7 


7670 


3780 





230 ± 


030 


0.199 ±0 


024 


0.08 


0.114 ± 0.030 


(3.80-) 


8 


6990 


3780 





231 ± 


029 


0.202 ±0 


023 


-0.21 


0.121 ± 0.031 


(3.9(7) 


9 


6900 


3750 





233 ±0 


029 


0.198 ±0 


024 


-0.18 


0.116 ±0.031 


(3.7cr) 


10 


7040 


3780 





234 ± 


029 


0.208 ±0 


023 


-0.15 


0.127 ± 0.031 


(4.1cr) 


(xl,7^) 


1 


6990 


3750 





199 ±0 


022 


0.160 ±0 


022 


-0.19 


0.080 ±0.031 


(2.6cr) 


2 


7770 


3750 





200 ± 


024 


0.160 ±0 


022 


-0.09 


0.077 ± 0.030 


(2.6cr) 


3 


6720 


3750 





203 ±0 


024 


0.166 ±0 


022 


-0.28 


0.088 ±0.031 


(2.80-) 


4 


6900 


3750 





204 ± 


024 


0.164 ± 


022 


-0.18 


0.083 ±0.031 


(2.7(t) 


5 


6760 


3750 





206 ±0 


027 


0.165 ±0 


022 


-0.22 


0.085 ±0.031 


(2.8(t) 


6 


6950 


3750 





206 ±0 


024 


0.166 ±0 


023 


-0.19 


0.085 ±0.031 


(2.70-) 


7 


6850 


3750 





206 ±0 


024 


0.166 ±0 


022 


-0.22 


0.085 ±0.031 


(2.8(t) 


8 


7930 


3750 





209 ±0 


027 


0.167 ± 


022 


-0.12 


0.085 ± 0.030 


(2.80-) 


9 


6590 


3750 





209 ±0 


028 


0.166 ±0 


023 


-0.23 


0.092 ± 0.032 


(2.9cr) 


10 


6990 


3780 





209 ±0 


027 


0.177 ± 


022 


-0.19 


0.097 ± 0.030 


(3.20-) 




1 


5580 


6330 





146 ±0 


021 


0.083 ±0 


022 


0.43 


-0.004 ± 0.018 


(0.2O-) 


2 


3980 


4140 





147 ± 


017 


0.076 ± 


025 


0.56 


-0.002 ±0.017 


(O.lo-) 


3 


6330 


5580 





149 ±0 


021 


0.086 ±0 


022 


0.46 


0.002 ±0.018 


(O.lcr) 


4 


5730 


6370 





149 ±0 


019 


0.087 ± 


022 


0.41 


0.002 ±0.019 


(O.lcr) 


5 


6370 


5580 





150 ± 


019 


0.094 ± 


021 


0.42 


0.009 ±0.018 


(0.5O-) 


6 


5580 


6290 





150 ± 


022 


0.076 ±0 


025 


0.39 


-0.012 ±0.021 


(0.6cr) 


7 


5730 


6330 





150 ± 


019 


0.087 ± 


023 


0.43 


-0.000 ±0.019 


(O.Ocr) 


8 


5580 


6420 





151 ± 


021 


0.093 ±0 


021 


0.43 


0.003 ±0.019 


(0.2O-) 


9 


5690 


6330 





151 ± 


021 


0.085 ±0 


022 


0.18 


0.009 ± 0.022 


(0.4O-) 


10 


6370 


5730 





151 ± 


018 


0.083 ±0 


024 


0.64 


-0.002 ±0.019 


(O.lcr) 


(a;i,c,7e'=) 


1 


3780 


4580 





130 ± 


017 


0.050 ±0 


029 


0.81 


-0.020 ±0.012 


(1.7cr) 


2 


3610 


4890 





132 ±0 


015 


0.056 ±0 


028 


0.71 


-0.016 ±0.014 


(1.2(t) 


3 


5360 


6900 





132 ±0 


017 


0.062 ±0 


025 


0.84 


-0.014 ± 0.010 


(1.4ct) 


4 


4760 


3800 





132 ±0 


020 


0.065 ±0 


025 


0.84 


-0.012 ±0.010 


(1.2ct) 


5 


7980 


8140 





133 ± 


015 


0.062 ±0 


025 


0.69 


-0.021 ± 0.015 


(1.4(t) 


6 


7720 


8200 





133 ±0 


015 


0.061 ± 


026 


0.94 


-0.013 ±0.008 


(I.60-) 


7 


6080 


6210 





133 ±0 


020 


0.054 ± 


029 


0.79 


-0.019 ±0.013 


(1.5cr) 


8 


3610 


4170 





133 ±0 


015 


0.064 ± 


025 


0.81 


-0.009 ±0.010 


(0.9(t) 


9 


4280 


3780 





133 ±0 


015 


0.062 ±0 


026 


0.80 


-0.011 ± 0.011 


(l.Ocr) 


10 


3610 


4230 





133 ±0 


017 


0.062 ±0 


024 


0.73 


-0.016 ±0.013 


(1.2(t) 


(xi,c) 








0.144 ± 0.019 


0.076 ± 0.023 









4.10.3 Model 2: 7^ and xi 

Combining a flux ratio with the SALT2 stretch parameter xi 
(Equation [5| also leads to no improvement over the {x\,c) 
model. The WRMS and o-prod values are smaller than when 
using the 7?.-only model, but they are still significantly larger 
than those from the standard {xi , c) model. We again point 
out that [Blondin et al.| ( |2011[ ) found that the (xi , TZ) model 
did better than the {xi, c) model. 

The best flux ratio using the {xi,TZ) model, 
7?. (6990/3750), is (like the 7?.-only model) not correlated 
with xi, but strongly correlated with c (coefficient of 0.83). 
The ratio TZ (6990/3750) is therefore a proxy for c, and so it 
is unsurprising that this is the top-ranked ratio in a model 
employing only a ratio and xi. Using this model [Blondin| 



et al. ( 2011 1 again found that the ratio 7^ (6630/4400) was 
best and, as mentioned above, it was similarly correlated 
with SALT2 c. Nearly all of the top ten ratios for the 7?.- 
only and {xi,TZ) models have very similar numerator and 
denominator wavelengths, and the difference in wavelength 
between the two fluxes is sig nificant (3000-4000 A). This 
again supports the idea presented above that the top-ranked 
flux ratios for these two models are effectively proxies for 
colour. 



4.10.4 Model 3: TZ" and c 

Some of the top-ranked flux ratios with the SALT2 colour 
parameter c (Equation m are consistent with the results 
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when using the (xi , c) model. This lack of improvement is 
once again at odds with what has been seen befo re ([Blondin 
et al. 20111. However, as mentioned in Section 4.10.2 the 



apparent improvement is likely due to the smaller sample 



size used by (Blondin et al. 20111 and the relatively poor 



performance of the {xi,c) model. 

Many ratios appear to have large correlations between 
the correction terms and uncorrected residuals, but the low- 
est WRMS values are tightly clustered in wavelength space. 
This is also apparent in Table [l] where all but one of the 
top ten ratios for the (c, TZ'^) model involve wavelengths near 
5600 A and 6300 A. Flux ratios which include these wave- 
lengths are effectively the same as 5R(SiS), the SiS ratio (Sec- 
tion ■ 



4.71 



The SiS ratio was shown above to be anticorrelated with 
Ami5(_B) (as well as correlated with xi). The best flux ratio 
using the (c, T?.'^) model, 7?. (5580/6330), is strongly corre- 
lated with xi (correlation coefficient of 0.83) and effectively 
uncorrelated with c (correlation coefficient of —0.18). There- 
fore, the top-ranked ratio for this model can be thought of 
as equivalent to 5R(SiS) and/or xi. Furthermore, the lack 
of correlation between TZ'' (5580/6330) and c implies that 
dereddening the data using the SALT2 c and colour law is 
working as intended. 

With the (c, TZ") model, [Blondin et"aL] ( |2011[ ) showed 
that their top-ranked ratio was TZ"^ (6420/5290). This is sim- 
ilar to the reciprocal of the best ratio found with the BSNIP 
data; thus, it is not surprising that they find a strong an- 
ticorrelation between 7?. (6420/5290) and xi. Their ratio 
also leads to a larger decrease in WRMS (~ 15 per cent. 



Blondin et al. 20111. Bailey et al. (20091 again showed an 



improvement when using their top ratio for this model, 
TZ'^ (6420/5190), which is also very close to the best ratio 
of Blondin et al. (20111, as well as the reciprocal of the 
best ratio presented here. Despite this similarity in wave- 
length space, the top-ranked ratios for this model from[BaF| 
ley et aLlpoM]) and|Blondin et al.|(|201l"|) yield WRMS val- 



ues of ~0.21 mag when using the BSNIP data. 



4.10.5 Model 4: TZ" , Xi, and c 

Quite significant improvements over the standard {xi,c) 
model are seen when using the top-ranked fiux ratios along 
with both xi and c (Equation [7|. In fact, the top ten ratios 
lead to improvements at about the l-2cr level (see Table [T]). 
The {xi, c, TZ") model and its best ratio, W (3780/4580), de- 
crease the WRMS by ~10 per cent and CTpred by ~34 per cent 
from the {xi,c) model. The Hubble diagram for the top- 
ranked fiux ratio using the {xi,c,TZ'^) model is shown in 



Figure 30 along with the WRMS for the best model (the 
grey band). 

The decrease in WRMS using this model is smaller 
than previously seen, but the decrease in Upred is larger and 
the overall improvement is at a higher significance; [Blondin] 
|et al.| ( |2011[ ) found a 20 per cent decrease in WRMS with 
a 1.6(7 significance. The fact that the current dataset has 
nearly three times as many objects as the sample studied 
by [Blondin et al!| ( |2011[ ) is likely the reason why the models 
using BSNIP spectra are showing less overall improvement 
over the standard {x\, c) model, but for the {xi, c, 7?.°) model 
the improvements are more significant. 

Most ratios have low values of WRMS and extremely 




R'(3780/4580) 



0.02 



0.04 



0.06 



0.08 



0.10 



Figure 30. Hubble diagram residuals versus ^CMB for the best 
flux ratio, 7^<= (3780/4580), using the (a•l,c,7^'=) model (Equa- 
tion [tJ. The grey band is the WRMS for the model. Colours and 
shapes of data points are the same as in Figure |5| 



large correlations between the correction terms and uncor- 
rected residuals, implying that the {xi,c,TZ'^) model is per- 
forming better overall than any of the other models investi- 



gated. This is consistent with Blondin et al. ( 2011 1, although 



they found no fiux ratio to have a strong correlation between 
the correction terms and the uncorrected Hubble residuals. 
The wavelengths of six of the top ten ratios for this model 
are near ~3750 A and ~4550 A with wavelength baselines of 
~800 A. These approximately correspond to the midpoint 
of the Ca II H&K feature and the border between the Mg II 
and Fe II complexes. 

Figure [sT] shows the best fiux ratio using the (a;i, c, VJ^) 
model, 7^'' (3780/4580), versus the SALT2 parameters xi 
and c. There is a slight correlation of this ratio with xi 
(correlation coefficient of 0.46) and effectively no correla- 
tion with c. Since this ratio is essentially uncorrelated with 
both SALT2 parameters and it decreases the WRMS at the 
~2cr level, it yields useful information about each SN beyond 
light-curve stretch and colour. It is intriguing that this new 
information is found at the blue end of the optical range, 
since there is evidence that spectral features in this region 
do not correlate with light-curve parameters, yet contain 
information related to SN la luminosity (e.g 
|2008| ) 



Foley et al. 



Note that Blondin et al. (20111 found that fiux ratios 
-5300 A and baselines of < 400 A 
, TZ'^) model. Using the 



with wavelengths near 
gave the best results for the (x\ , c, 
best ratio from their study, TZ" (5690/5360), the BSNIP data 
yield a WRMS of ~0.15 mag, which is not as good as our 
top ten ratios. 



5 CONCLUSIONS 

This is the third paper in the BSNIP series and presents 
a comparison between spectral feature measurements and 
photometric properties of 108 low-redshift {z < 0.1) SNe la 
within 5 d of maximum brightness. The spectral data all 
come from BSNIP I, and the photometric data come mainly 
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Figure 31. The top-ranked flux ratio, 11" (3780/4580), using the 
(xi,c, T?.'^) model (Equation [tJ versus the SALT2 parameters xi 
and c. Colours and shapes of data points are the same as in Fig- 
ure [5] 



from the LOSS sample and are published by lGaneshalingam] 
|et al.| p010[ ). The details of the spectral measurements can 
be found in BSNIP II, and the light-curve fits and photo- 
metric parameters are in Ganeshalingam et al. (in prepara- 
tion). A combination of light-curve parameters (specifically 
the SALT2 stretch and colour parameters X\ and c) and 
spectral measurements are used to calculate distances to 
SNe la. We then compare the residuals from these models to 
the standard model which only uses light-curve stretch and 
colour. Future BSNIP papers will incorporate host-galaxy 
properties and SN spectra at later epochs into the analysis 
presented here. 



5.1 Summary of Investigated Correlations 

The velocity gradient ( [Benetti et al.|[2005| ) is compared to 
the light-curve width, and it is shown that, as in BSNIP II, 
the various classifications based on the value of v overlap sig- 
nificantly. Similarly, velocities at maximum brightness (uq) 
are compared to photometric observables and classifications 



based on velocity gradient, and there is a large amount of 
overlap in all of these parameters as well. In earlier work, HV 
and HVG objects have been used almost interchangeably, as 
have normal velocity and LVG objects (e.g., Hachinger et al 



2006[ [Pignata et al.||2008| |Wang et al.||2009[ ). However, the 



analyses of BSNIP II and this work show that these associ- 
ations are not as distinct as previously thought. 

The measured velocities of the Si II A6355 and Si II 
A5972 features near maximum brightness are uncorrelated 
with observed {B — V)^^^. Furthermore, the HV objects 
and normal-velocity objects have similar distrubutions of 
observed [B — V)^^^ values, though the HV objects tend to 
have slightly larger observed (B — ^)jnax- When distances 
to SNe la are computed using light-curve width {xi) and 
colour (c) parameters, no significant improvement in the ac- 
curacy of the distances is found when the velocity of the Si II 
A6355 feature is added. This is opposite to what was seen 



by Blondin et al. (2011 1. Furthermore, models involving x\, 
c, and the velocity of the Ca II H&K feature fare as poorly. 
However, the velocity of the S II "W," when used in con- 
junction with x\ and c, leads to a decrease in the WRMS of 
the distances at the 1.8cr level. Despite what was seen in the 



analysis of Blondin et al. ( 2011 1, the use of relative depths of 



none of the features analysed here leads to an improvement 
in the Hubble residuals. 

The pEW of the Si II A4000 feature is strongly anti- 
correlated with x\, which confirms many previous studies 



(Arsenijevic et al. 2008 Walker et al. 2011 Blondin et al 



2011 Nordin et al. 2011b Chotard et al. 20111. Further 



more, when using a model that includes x\, c, and the pEW 
of this feature, the residuals are as low as when using the 
standard (xi , c) model. The pEWs of the Mg II and Fe II 
complexes are both correlated with c, yet when using either 
of these (as a proxy for c) combined with the pEW of the 
Si II A4000 feature (as a proxy for 2:1), the Hubble diagram 
residuals are significantly larger than when simply using xi 
and c. 

The pEWs of both Si II A5972 and Si II A6355 are well 
correlated with xi and correlated with c, but the use of the 
Si II A5972 pEW does not improve distance calculations. 
However, using the Si II A6355 pEW (along with xi and c) 
leads to an improvement in the WRMS residuals at the 1.2a 
level. Finally, the Ca II near-IR triplet is correlated with c 
and MLCS2k2 light-curve width parameter A. This feature 
and the O I triplet have not been investigated thoroughly in 
studies similar to this one since other large SN la spectral 
datasets often do not include these spectral regions. 

The Si II ratio, used as a luminosity indicator previously 



(e.g., Nugent et al.|1995 Benetti et al.|2065 Hachinger et al 
|2006[ ), is found to be well correlated with Amis. However, we 
caution that at a given value of Amis , there can exist various 
spectroscopically classified subtypes of SNe la. This ratio is 
also found to be correlated with observed (B — V) , which 
has been seen in other work ( Altavilla et al.|2009 |. We also 
show that the Si II ratio is not an accurate proxy for xi when 
calculating distance moduli. A model using c and 5R(Si II) 
performs significantly worse than the usual {xi,c) model, 
contrary to the conclusion of [Blondin et al. (2011 1. 

On the other hand, the Ca II ratio is found to be a 
good indicator of light-curve width, as it is well correlated 
with the MLCS2k2 A parameter as well as with Amis. The 
BSNIP data also indicate that the SiS ratio is correlated 



28 Silverman, et al. 



with both xi and c and distance models using 5R(SiS) with 
just c or with both c and xi perform as well as the standard 
{xi,c) model. Finally, we confirm the results of iHachinger, 



et al. ( 2006 \ that the SSi and SiFe ratios are both accurate 



luminosity indicators, as they are both well correlated with 

Amis- 

Following Bailey ct al.' ('2009^ and [Blondin et aT] poTTj) , 



we calculate Hubble diagram residuals using models which 
include combinations of the usual light-curve parameters 
(width and colour) and arbitrary sets of flux ratios. A total 
of 17,822 different ratios of fluxes are used alone, with xi, 
with c, and with both xi and c to investigate whether any 
of these models might improve the accuracy of SN la dis- 
tance measurements. No models utilising only a flux ratio or 
a flux ratio and xi are found to decrease the Hubble residu- 
als. A handful of models using a flux ratio and c are seen to 
perform as well as the standard {xi , c) model. Interestingly, 
most of these best ratios are extremely close to the SiS ratio 
mentioned above. 

These results differ from those of the previous studies 
of [Bailey et~aL] ( |2009| and [Blondin et al.| ( |20TT| ), both of 
which found that flux ratios alone or in conjunction with 
light-curve information would usually perform better than 
the (xi, c) model. This is likely due to the fact that the num- 
ber of objects used in the current analysis is 2-3 times larger 
than what was used in these past studies and our "standard 
model" (without a spectral indicator) already performs sig- 



niflcantly better than that of Bailey et al. ( 2009 \ or Blondin 



et al. (20111. The differences may also be caused by the 



larger number of spectroscopically peculiar SNe la in the 
BSNIP sample. 

Finally, when combining a flux ratio with both xi and 
c, our top-performing ratio, 7?. (3780/4580), decreases the 
Hubble residuals by 10 per cent, which is significant at the 
2(7 level. The WRMS of the residuals using this model is 
0.130 ± 0.017 mag, as compared to 0.144 ± 0.019 mag when 
using the same sample with the standard {xi,c) model. 
This Hubble diagram has one of the smallest scatters ever 
published and at the highest significance ever seen in such 
a study. The wavelengths involved in most of the best- 
performing ratios in the {xi,c,TZ'^) model approximately 
correspond to the midpoint of the Ca II H&K feature and 
the border between the Mg II and Fe II complexes. This sup- 
ports previous work which has shown that near-UV spectra 
of SNe la contain information related to SN la luminosity 
which is not necessarily captured in the photometry (e.g., 
Foley et al.|2008| ). 



5.2 The Future? 

New large-scale surveys are already obtaining SN la data, 
and they are observing to higher redshifts and gathering 
larger amounts of data than what is in the BSNIP sample 



(e.g., [Rau et al.[[2009| [Law et al.[[2009[ [Kaiser et aL][2002 1 
Even larger surveys at even higher redshifts are also planned 
(e.g., LSST, WFIRST). Many more SNe la will be discovered 
than can be rigorously observed; we are quickly entering the 
age of SN research where we are limited by the follow-up 
observations. Thus, there must be significant effort put forth 
to determine the most efficient way to monitor and utilise 
such vast quantities of objects. That is one of the major 
goals of BSNIR 



Soon, for the vast majority of objects, there will only be 
(at best) a handful of photometric observations near maxi- 
mum brightness. Those, combined with a relatively low S/N 
spectrum near maximum, will likely be all the follow-up ob- 
servations we get. From the work presented here (and in 
BSNIP II) we have shown that there is still hope. 

The pEW of the Si II A4000 feature is a good indica- 
tor of light-curve width, and the pEWs of the Mg II and 
Fe II complexes are relatively good proxies for colour. Un- 
fortunately, the correlations between these spectral mea- 
surements and the corresponding photometric properties 
is not perfect, and distance calculations that employ only 
these spectroscopic measurements do not perform as well 
as the standard model which uses light-curve width and 
colour. However, this is still a promising avenue for fur- 
ther investigation using new datasets that are even larger 
than BSNIP. Other correlations that appear marginal in the 
BSNIP dataset, or models tested here that performed only 
equally as well as the usual {xi,c) model, should also be 
reexamined in the future. 

Occasionally, one will be fortunate enough to have suf- 
ficient photometric observations to produce a light curve 
for which SALT2 (or another light-curve fitter) is able to 
determine a width and colour. In these cases it appears 
that the light-curve parameters can be combined with a 
flux ratio from a spectrum near maximum brightness to 
improve the accuracy of SN la distances. The best ratios 
for this, as determined from the BSNIP data, are all near 
7^'=(~3750/~4550). 

This is all somewhat heartening for surveys that will 
discover and monitor SNe la at higher redshifts. Si II A4000, 
the Mg II and Fe II complexes, and 7^'=(~3750/~4550) aU 
involve spectral features which are toward the blue end of 
the optical range. This is critical for higher-z surveys since, 
as pointed out in BSNIP II, the red wing of the typical, 
near-maximum Si II A6355 feature becomes redshifted be- 
yond ~l/im for z > 0.6. Furthermore, as discussed multi- 
ple times in BSNIP II, measuring fluxes and pEWs directly 
from a spectrum is much easier and less reliant on smooth- 
ing models or functional form assumptions than velocities, 
for example. 

To quote the concluding paragraph of [Blondin et al.| 
(20111, "Do spectra improve distance measurements of 
SN la? Yes, but not as much as we had hoped." "Good" 
scientists are not supposed to let personal feelings or hopes 
interfere with their work. However, we have to agree with the 
authors of this quote both on the objective part (i.e., spec- 
tra do improve distance measurements), as well as on the 
subjective part (i.e., we hoped they would improve things 
even more). 
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